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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 

The  output  wavelength  of  PbSe  diode  lasers  operating  at  77°K  has  been  tuned  in  the  range 
from  7.  5  to  over  22 by  applying  hydrostatic  pressures  of  up  to  14  kbar.  Extension  to  wave¬ 
lengths  above  13m  was  made  possible  by  the  development  of  polycrystalline  germanium  and 
silicon  windows  capable  of  withstanding  the  high  pressures  involved.  These  windows  have 
been  tested  to  20  kbar  with  no  evidence  of  deterioration.  Germanium  windows  were  used  in 
the  wavelength  range  to  20m  in  order  to  eliminate  transmission  difficulties  due  to  lattice  ab¬ 
sorption  bands  in  silicon;  above  20 m>  where  germanium  has  lattice  absorption  bands,  silicon 
windows  were  used.  Our  range  of  tuning  has  been  limited  because  the  helium  pressure- 
transmitting  fluid  freezes  at  14.15  ±0.03  kbar  at  77°K.  Since  hydrostatic  pressure  cannot 
be  transmitted  reliably  through  solid  helium,  we  plan  to  extend  the  diode  emission  to  longer 
wavelengths  by  first  pressurizing  the  system  to  above  14  kbar  at  room  temperature  (where 
helium  is  still  a  gas),  and  then  cooling  it  to  77°K  and  below.  With  constant  current  through 
the  diode,  an  increasing  pressure  has  been  observed  to  produce  no  significant  change  in  the 
intensity  of  either  the  coherent  or  spontaneous  emission,  until  the  freezing  point  of  helium 
is  reached.  The  large  reduction  in  the  intensity  of  emission  that  occurred  when  the  helium 
froze  is  attributed  to  diffuse  scattering  of  the  radiation  by  the  solid  helium. 

A  new  technique  has  been  developed  for  growing  a  thin  layer  of  InSb  on  an  InSb  substrate  from 

a  saturated  In-InSb  melt.  This  type  of  growth  from  a  solution  has  been  previously  used  in 

the  fabrication  of  GaAs  and  Ge  p-n  junctions;  however,  at  the  lower  temperatures  required 

for  InSb  (about  300°C),  certain  modifications  were  made,  such  as  the  use  of  hot  liquid  stearic 

acid  which  covers  the  InSb  substrate  and  the  In-InSb  melt  during  the  growth  process.  By 

suitably  doping  the  melt,  uniform  planar  InSb  n+-p  junctions  were  grown  on  substrate  areas 
2 

of  about  2  cm  .  Such  InSb  junctions  have  been  used  in  large-volume  injection  lasers  which 
emit  coherent  light  in  a  direction  normal  to  the  junction  plane. 

Injection  of  carriers  into  InSb  n+pp+  structures  similar  to  those  used  as  bulk  lasers  has  been 

studied  by  determining  the  distribution  of  potential  and  electric  field  along  the  p-region  at 

various  current  levels.  Measurements  were  made  using  a  movable  potential  probe  with  the 

2 

sample  at  77°K.  These  measurements  indicate  that  at  currents  of  less  than  1  A/cm  the 
dominating  mechanism  is  the  injection  of  electrons  into  the  p-region  from  the  n+  contact,  which 
causes  a  saturation  of  electron  traps  and  an  increased  lifetime  in  the  p-type  bulk.  At  higher 
currents  an  appreciable  number  of  holes  is  injected  from  the  p+  region  in  order  to  neutral¬ 
ize  the  large  number  of  injected  electrons.  This  changeover  from  "single"  to  "double"  in¬ 
jection  produces  changes  in  the  potential  distribution  as  well  as  changes  in  the  region  of  bulk 
InSb  from  which  the  laser  radiation  is  emitted. 

Microwave  amplification  has  been  observed  in  bulk  n-GaAs  at  room  temperature.  The  sam- 

2 

pies  were  50-M-thick  rectangular  parallelepipeds  with  a  mobility  of  about  5000  cm  /V-sec 
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and  a  resistivity  of  approximately  lOOohm-cm.  The  frequency  at  which  the  maximum  gain 
occurs  agrees  fairly  well  with  the  frequency  at  which  Gunn  oscillations  occur  in  lower  resis¬ 
tivity  samples  (~1  ohm-cm)of  the  same  thickness,  but  the  average  applied  electric  field  nec¬ 
essary  for  gain  is  about  7000  V/cm  —  about  twice  the  value  usually  required  for  Gunn  oscil¬ 
lations.  The  amplifier  sample  exhibits  only  a  slight  tendency  toward  current  saturation  with 
increasing  voltage,  whereas  the  oscillator  sample  begins  to  saturate  even  below  the  thresh¬ 
old  of  oscillation,  and  above  threshold  exhibits  a  range  of  current  oscillation  relatively  in¬ 
dependent  of  voltage.  The  existence  of  a  positive  resistance  in  the  amplifier  samples  for 
average  applied  electric  fields  much  higher  than  the  field  at  which  Gunn  oscillation  normally 
occurs  suggests  that  the  usual  domain  formation  does  not  take  place  in  these  high-resistivity 
samples.  This  agrees  with  recent  calculations  of  McCumber  and  Chynoweth  who  have  found 
that,  for  sufficiently  small  products  of  carrier  concentration  and  sample  length,  the  trans¬ 
ferred  electron  effect  leads  to  a  negative  conductance  at  frequencies  near  the  reciprocal  of 
the  transit  time,  even  though  the  DC  solution  displays  only  positive  resistance. 

II.  LASER  RESEARCH 

Absolute  (nonconvective)  instabilities  were  not  found  from  an  analysis  of  the  dispersion  re¬ 
lation  for  Stokes  and  anti-Stokes  modes  in  an  infinite  Raman-active  medium.  Absolute  in¬ 
stabilities  were  found  at  zero  frequency  and  second  harmonic  frequency  of  the  laser,  but  in¬ 
accuracy  of  the  dispersion  relation  at  these  frequencies  makes  dubious  the  physical  realiza¬ 
tion  of  these  instabilities,  which  have  approximately  millisecond  growth  times. 

Calculations  have  been  made  of  the  focusing  distance  and  intensification  of  a  self-focused 
high-power  laser  beam.  Self-focusing  is  due  to  an  increase  of  the  index  of  refraction  of  the 
propagation  medium  which  is  proportional  to  the  local  laser  intensity.  Computer  solutions 
were  obtained  from  an  approximate  nonlinear  wave  equation,  which  is  valid  when  the  inten¬ 
sity  varies  slowly  over  a  wavelength,  for  an  input  beam  with  a  plane  wavefront  and  a  Gauss¬ 
ian  radial  intensity  profile.  Typically,  a  1-MW  beam  of  2-mm  diameter  would  focus  itself  in 
about  50  cm  of  CS2* 

Dielectric  breakdown  in  water  produced  by  a  focused  Q-switched  ruby  laser  resulted  in  an 

o 

unexplained  blue  emission  with  wavelengths  ranging  from  3200  to  4500  A,  with  a  minimum 

o 

at  3600 A.  The  emission  lasted  20 nsec  or  less,  aside  from  a  long  but  weak  afterglow. 

Fluorescence  from  limonene,  possibly  produced  by  multiple-photon  absorption  from  a  plane- 
polarized  ruby  laser  beam,  showed  little  depolarization  due  to  Brownian  rotation  of  the  ex¬ 
cited  molecules  during  their  radiative  lifetime.  Considering  the  value  of  the  rotation  relaxa¬ 
tion  time,  one  estimates  that  the  radiative  lifetime  is  less  than  0.1  nsec,  which  is  consistent 
with  direct  measurements  of  the  fluorescent  decay. 

Analysis  of  the  symmetry  of  molecular  polarizability  tensors  shows  that  no  incoherent  emis¬ 
sion  of  sum  frequency  light  will  arise  from  density  fluctuations  unless  the  medium  is  opti¬ 
cally  active  and  unless  the  two  primary  frequencies  are  different;  otherwise,  the  sum  fre¬ 
quency  emission  may  arise  from  orientation  fluctuations. 
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Stimulated  Stokes  and  anti-Stokes  lines  have  been  obtained  in  quartz,  using  a  high-power 
ruby  laser  as  a  source.  As  part  of  an  investigation  of  the  interference  between  Brillouin  and 
Raman  effects  in  solids,  details  of  the  spectrum  have  been  compared  with  spontaneous  Raman 
scattering  in  quartz. 

The  Bose-Einstein  distribution  has  been  confirmed  for  the  photoelectron  counts  from  a  photo¬ 
multiplier  which  is  illuminated  by  a  gas  laser  operating  just  below  threshold,  as  expected 
for  a  narrow-band  Gaussian  source  when  the  counting  interval  is  less  than  the  inverse  band¬ 
width  of  the  light. 

III.  MATERIALS  RESEARCH 

Liquidus  temperatures  for  the  portions  of  the  Hg-Te  system  from  3  to  25  and  from  50  to  100 
atomic-percent  tellurium  have  been  determined  by  differential  thermal  analysis.  These  re¬ 
sults  have  been  combined  with  liquidus  points  obtained  earlier  from  partial  pressure  data  to 
establish  the  Hg-Te  phase  diagram. 

The  pressures  of  mercury  and  selenium  in  the  vapor  over  Hg-saturated,  Se-saturated,  and 
congruently  subliming  HgSe(c)  between  450°  and  800°C  have  been  determined  by  measuring 
the  optical  absorption  of  the  vapor.  The  pressure  data  have  been  used  to  obtain  the  first  re¬ 
ported  values  for  the  standard  Gibbs  free  energy  of  formation  of  HgSe(c). 

At  sufficiently  low  temperatures,  the  high-pressure  phase  of  InTe  can  be  retained  at  atmos¬ 
pheric  pressure  as  a  metastable  phase.  The  heat  evolved  when  this  metastable  phase  is 
transformed  at  atmospheric  pressure  into  the  stable  low-pressure  form  has  been  measured 
by  both  metal-solution  and  differential -scanning  calorimetry.  The  heat  of  transformation  is 
0.44  kcal/g-atom. 

High-pressure  phases  can  also  be  retained  at  atmospheric  pressure  in  InSb(II)  _x  (2  In)x 
samples  kept  at  sufficiently  low  temperatures.  Qualitative  results  on  the  amount  of  heat 
liberated  when  such  samples  are  transformed  at  atmospheric  pressure  indicate  the  existence 
of  a  high-pressure  phase  different  in  structure  from  either  InSb(II)  or  indium.  This  phase 
lies  between  x  =  0.  17  and  x  =  0.  60,  but  its  exact  composition  has  not  been  determined. 

The  temperature  at  which  the  phase  of  MnAs  with  B8  structure  is  transformed  into  the  high- 
temperature  B31  form  has  been  measured  as  a  function  of  hydrostatic  pressure.  The  tran¬ 
sition  temperature  decreases  from  44°C  at  atmospheric  pressure  to  — 72°C  at  4.5kbar.  The 
functional  dependence  of  transition  temperature  on  pressure  indicates  that  the  B3I  phase  is 
stable  at  all  temperatures  for  pressures  exceeding  4.5kbar. 

Electrical  measurements  have  been  made  on  T^O^  single  crystals  containing  considerably 
less  nitrogen,  the  principal  impurity,  than  earlier  samples.  In  contrast  to  previous  results, 
both  resistivity  and  Hall  coefficient  depend  strongly  on  orientation.  The  magnetoresistance, 
A P/PQ.  is  2.5  at  4.2°K  and  170  kG,  compared  with  the  previous  maximum  value  of  0.12. 

The  concept  of  a  Jahn-Teller  theorem  for  narrow-band  electrons  has  been  formalized  and 
applied  to  the  interpretation  of  the  localized  vs  collective  character  of  outer  d-electrons  in  a 
wide  range  of  transition-metal  compounds. 
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Wet  chemical  methods  have  been  developed  for  determining  cobalt  and  the  rare  earth  element 
in  the  lanthanum  and  yttrium  cobalt  oxides,  and  for  determining  all  three  major  components 
in  iron-nickel-copper  thin  films. 


IV.  PHYSICS  OF  SOLIDS 

Reflectivity  measurements  in  ReO^,  which  were  previously  reported  in  the  range  from  0.2 
to3.8eV,  have  now  been  extended  to  12 eV.  Results  indicate  a  Drude-like  behavior  in  the 
energy  region  below  2eV  and  peaks  at  4.6  and  9.6  eV  arising  from  interband  transitions.  A 
Kramers-Kronig  analysis  of  the  reflectivity  has  been  used  to  obtain  the  real  and  imaginary 
parts  of  the  dielectric  constant,  and  a  further  separation  of  the  dielectric  constant  compo¬ 
nents  into  free-electron  and  interband  contributions  has  been  made. 

Concurrently,  a  tight  binding  scheme  is  being  used  to  investigate  the  electronic  band  struc¬ 
ture  of  ReO^.  The  energy  levels  have  been  found  to  separate  into  cr  and  7r  bonding  and  anti¬ 
bonding  bands,  and  into  a  set  of  nonbonding  n  bands.  Estimates  of  the  various  parameters 
are  being  made,  using  the  augmented  plane-wave  method  to  calculate  the  energy  values  at 
the  high  symmetry  points  of  the  Brillouin  zone.  This  calculation  will  be  used  to  predict  opti¬ 
cal,  magnetic,  and  transport  properties. 

Experimental  investigation  of  the  band  structure  of  GaSe  is  continuing.  A  reduced  effective 
mass  m*  =0.l4m  and  an  energy  gap  (b^  =  2. 130  eV  have  been  obtained  from  more  highly 
resolved  oscillatory  magnetoabsorption  data  at  1.5°K  in  the  Faraday  configuration.  The  data 
in  the  Voigt  configuration  are  smaller  in  amplitude  and  more  complex,  suggesting  that  more 
than  one  transition  is  involved.  In  the  zero-field  data,  three  clearly  defined  exciton  peaks 
have  been  resolved;  these  have  been  used  to  calculate  a  binding  energy  of  21  meV  and  a  re¬ 
duced  effective  mass  m*  =  0.17  m,  which  is  consistent  with  the  results  of  the  magneto-optical 
data.  Several  models  for  GaSe  which  are  consistent  with  the  experimental  data  have  been 
proposed. 

GaSe  is  also  being  examined  in  the  microwave  region  at  70  GHz  in  order  to  observe  intraband 
cyclotron  resonance.  So  far,  a  well-resolved  magnetoabsorption  peak  corresponding  to  an 
effective  mass  m*  =  0.02  m  and  an  o>t  ~  1.5  has  been  observed,  but  only  in  one  sample. 

Another  microwave  resonance  experiment,  on  p-type  PbSe,  is  continuing.  Marked  structure 
has  been  obtained  below  3kG  and  identified  with  dielectric  anomalies;  this  behavior  is  con¬ 
sistent  with  that  predicted  on  the  basis  of  a  classical  skin  effect  theory,  using  an  infinite 
relaxation  time  and  values  of  effective  mass  and  anisotropy  obtained  from  previous  Shubnikov- 
de  Haas  measurements  on  the  same  sample. 

A  new  approach  to  the  tight-binding  technique  for  phenomenologically  obtaining  the  energy 
bands  in  solids  has  been  developed.  This  method  involves  a  Fourier  expansionfor  the  energy 
bands  in  k-space  in  which  the  symmetry  of  the  lattice  and  the  degeneracies  of  the  bands  at 
high  symmetry  points  are  included.  The  Fourier  coefficients  of  the  expansion  are  evaluated 
from  experimentally  observed  parameters.  This  approach  is  being  tested  on  the  metal, 
aluminum,  and  the  semimetal,  antimony. 


An  electron -beam  technique  has  been  developed  for  the  vapor  deposition  of  piezoelectric 
semi- insulating  CdS  thin  films  on  various  substrates.  By  proper  orientation  of  the  electric 
field,  either  longitudinal  or  transverse  hypersonic  waves  can  be  excited.  Acoustically  active 
films  have  been  deposited  on  metals,  semiconductors,  and  insulators;  pulse-echo  techniques 
have  been  used  to  test  the  acoustical  activity  of  a  film.  The  range  of  frequency  excitation  in 
the  longitudinal  mode  for  active  films  varied  from  12.4  MHz  for  thick  films  deposited  on 
aluminum  to  70  GHz  for  thin  films  deposited  on  a  z-cut  quartz  rod. 

The  study  of  the  magnetic  properties  of  the  spinels  AC^X^,  where  A  is  nonmagnetic  zinc 
or  cadmium  and  X  is  oxygen,  sulfur,  or  selenium,  is  continuing.  Magnetic  measurements 
on  CdC^O^,  an  antiferromagnet  with  a  Neel  temperature  T^  =  9°K,  have  been  completed; 
CdC^S^  has  been  found  to  be  ferromagnetic  with  a  Curie  temperature  T^  =  86 °K.  It  has 
been  found  that,  in  addition  to  the  large  dependence  on  lattice  size  for  a  given  type  of  A-site 
ion,  T^,  the  asymptotic  Curie  temperature,  in  this  entire  series  of  compounds  depends 
strongly  on  whether  the  A-site  ion  is  zinc  or  cadmium.  Since  T^  is  proportional  to  the  sum 
of  all  magnetic  exchange  interactions,  this  indicates  that  the  diamagnetic  A-site  cation  signifi¬ 
cantly  affects  these  interactions.  Specifically,  substituting  zinc  for  cadmium  tends  to  de¬ 
crease  the  ferromagnetic  nearest-neighbor  interaction,  or  to  increase  the  antiferromagnetic 
next -nearest-neighbor  interaction,  or  both. 


In  order  to  resolve  inconsistencies  in  MnC^O^  between  the  nuclear  magnetic  resonance  and 
neutron  diffraction  results  for  the  manganese  cone  angle,  these  properties  were  reinvesti¬ 
gated.  The  new  data  are  in  general  agreement  with  the  earlier  findings,  but  indicate  that 
the  Mn4^  ions  possess  a  reduced  magnetic  moment  of  about  4.3  rather  than  the  previ¬ 
ously  assumed  value  of  5.0^.  However,  while  use  of  this  reduced  moment  improves  the 
agreement  between  the  theoretical  and  experimental  neutron  diffraction  patterns,  the  dis¬ 
crepancy  with  the  nuclear  magnetic  resonance  result  remains. 


An  experimental  study  of  the  temperature  variation  of  the  spin  wave  dispersion  relation  e,  = 

2  4  °  * 

Dk  —  Fk  has  indicated  that  for  a  63%Ni-37%  Fe  film  6950  A  thick,  both  the  quartic  and  the 

quadratic  coefficients  vary  as  T  '  up  to  80°K  and  approximately  as  T  '  up  to  298 °K.  A 
theoretical  justification  for  the  temperature  dependence  of  D  and  F  in  the  spin  wave  disper¬ 
sion  relation  can  be  simply  made  by  extending  a  model  proposed  by  W.  Marshall  for  treating 

a  system  of  interacting  spins  which  gave  the  dependence  of  D.  Higher-order  terms 

5/2 

occurring  in  D,  as  well  as  the  temperature  dependence  of  F,  predict  a  T  '  dependence  for 
7/2 

F,  a  higher-order  T  '  term  for  D,  and  a  relation  between  the  three  temperature  coeffi¬ 
cients  of  D  and  F. 


An  expression  has  been  derived  that  relates  the  fluctuation  scattering  of  light  expected  from 
a  ferromagnet  to  the  observed  Faraday-  or  Kerr-effect  coefficient  and  to  the  spin-spin  corre¬ 
lation  function.  An  estimate  is  made  of  the  fluctuation  scattering  to  be  anticipated  near  the 
Curie  temperature. 

Experimental  methods  of  observing  zero  sound  modes  in  metals  are  considered.  It  is  con¬ 
cluded  that  a  mode  characterized  by  a  paramagnetic  spin  wave  is  likeliest  to  exist,  and  should 
most  readily  be  observed  in  metals  which  show  superconductivity.  It  is  argued  that  these 
zero  sound  modes  may  exist  in  the  superconducting  state  as  well. 


viii 


A  simple  and  general  argument  has  been  found  for  considering  pseudomomentum  for  an  en¬ 
semble  of  interacting  phonons.  It  is  shown  that  the  anticipated  selection  rule  on  pseudo¬ 
momentum  holds  for  any  process  produced  by  phonon-phonon  interaction. 

The  energy  spectrum  of  a  polaron  has  been  considered  in  the  region  where  its  excitation 
energy  is  close  to  that  of  an  optical  phonon.  Both  a  variational  and  weak  coupling  calculation, 
which  are  consistent  with  each  other,  have  been  performed.  It  is  found  that  the  dispersion 
relation  in  the  weak  coupling  theory  becomes  flat  as  the  polaron  momentum  approaches  the 
value  at  which  the  polaron  excitation  energy  and  the  optical  phonon  energy  become  equal. 

A  treatment  has  been  formulated  of  the  kinetic  equation  for  a  system  of  electrons  under  the 
action  of  arbitrary  space-  and  time-varying  fields  and  a  random  distribution  of  impurities. 
The  results  can  be  used  to  derive  a  Boltzmann  equation  for  the  system.  Kinetic  equations 
for  other  situations,  including  the  presence  of  a  quantizing  magnetic  field,  can  be  readily 
obtained. 

The  interaction  of  photons  with  charged  particles  in  a  plasma  modifies  the  dispersion  re¬ 
lation  for  electromagnetic  propagation  in  a  well-known  way.  It  is  shown  that  this  result  can 
be  obtained  quite  simply  from  quantized  field  theory,  and  that  it  arises  because  of  the  A 
term  in  the  electromagnetic  energy,  averaged  over  the  charged  particles  of  the  medium. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  EXTENSION  OF  PRESSURE  TUNING  OF  PbSe  DIODE  LASER  TO  22  m  * 

The  measurements  reported  in  the  previous  Solid  State  Research  Report  (1965:2)  have  been 
extended  to  higher  pressures  with  the  development  of  windows  made  of  high-resistivity  polycrys¬ 
talline  germanium  or  silicon  formed  in  the  shape  of  cylinders  12  mm  diameter  by  12  mm  long. 

The  ends  of  the  cylinders  were  optically  polished.  One  end  was  placed  over  a  3 -mm  hole  in  the 
pressure  vessel.  Under  these  conditions  all  windows  tested  held  to  20  kbar  with  no  evidence  of 
deterioration.  Germanium  windows  were  used  for  wavelengths  to  20  1 1  in  order  to  eliminate 
transmission  difficulties  due  to  lattice  absorption  bands  in  silicon  between  8  and  20  \±;  silicon 
windows  were  used  for  wavelengths  above  20  p,  where  germanium  has  lattice  absorption  bands. 

Development  of  these  windows  has  made  possible  an  extension  of  the  experiments  to  the 
freezing  pressure  of  helium  at  77°K,  which  was  found  to  be  14.15  ±  0.03  kbar.  The  correspond¬ 
ing  diode  emission  wavelength  was  22  |j..  Pressure  tuning  of  the  PbSe  diode  emission  wavelength 
from  7.5  to  22  p  is  shown  in  Fig.  1-1.  As  the  helium  freezes,  the  rate  of  change  in  the  wavelength 
of  the  emission  with  applied  pressure  decreases  abruptly,  making  possible  an  accurate  determi¬ 
nation  of  the  freezing  pressure  of  helium.  Since  the  pressure  cannot  be  transmitted  reliably 
through  the  solid  helium,  we  plan  to  extend  the  diode  emission  to  longer  wavelengths  by  pressuriz¬ 
ing  the  system  to  above  14 kbar  at  room  temperature  (where  helium  is  still  a  gas),  and  then  cooling 
it  to  77  °K  and  below.  Also,  in  our  preliminary  experiment,  once  the  helium  froze,  the  diode  laser 
emission  decreased  to  the  point  where  it  was  barely  discernible.  Since  tfre  emission  intensity 
returned  to  its  original  value  each  time  the  pressure  was  reduced  below  the  helium  freezing 
pressure,  the  diode  suffered  no  permanent  physical  damage  in  the  process.  The  emission  in¬ 
tensity  decrease  is  believed  to  be  due  to  scattering  of  the  radiation  by  the  solid  helium.  To 
minimize  scattering,  the  diode  package  has  been  modified  so  that  the  diode  emitting  edge  can  be 
placed  within  0.05  mm  of  the.  window  surface. 

J.  M.  Besson* 

A.  R.  Calawa 
W.  Pault 

B.  SOLUTION  REGROWTH  OF  PLANAR  InSb  LASER  STRUCTURES 
1.  Introduction 

The  growth  of  thin  crystal  layers  from  a  melt  was  first  used  by  Nelson*  as  a  technique  of 
producing  remarkably  flat  and  abrupt  p-n  junctions  for  germanium  tunnel  diodes  and  GaAs 
junction  lasers.  In  the  case  of  InSb,  this  process  is  of  particular  interest  for  growing  an  n-type 
layer  on  a  p-type  substrate,  because  of  the  difficulty  of  producing  such  layers  by  diffusion. 
Diffusion  of  n-type  impurities  such  as  selenium,  tellurium,  or  sulfur  is  impractical  because 
these  elements  tend  to  form  glass -like  compounds  on  the  substrate  surface,  and  their  diffusion 
coefficients  in  InSb  are  small  at  temperatures  below  the  melting  point  of  InSb  (525 °C).  By  the 


*  The  work  described  in  this  section  was  supported  in  part  by  the  AF  Avionics  Laboratory  Director’s  Fund,  Item  Nr.  65-86. 
t  Harvard  University. 
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Section  I 


Fig.  1-1 .  Variation  with  pressure  of  dominant  emission  modes  from  PbSe  diode  laser  at  77° K. 

solution-regrowth  technique,  an  n-type  layer  can  be  grown  easily  on  a  p-type  substrate  by  simply 
adding  an  n-type  dopant,  such  as  tellurium,  to  the  melt. 

2.  Method  and  Results 

As  shown  in  Figs.  I-2(a-b),  the  apparatus  consists  of  a  boron-nitride  container  which  is 
heated  from  the  bottom  by  passing  current  through  a  graphite  heater  strip  and  which  can  be 
tilted  at  various  times  during  a  run.  One  side  of  the  boat  contains  the  InSb  substrate,  while 
the  opposite  side  contains  indium  and  a  sufficient  amount  of  InSb  to  produce  a  nearly  saturated 
In-InSb  solution  at  the  maximum  temperature  of  about  300° C.  When  n-type  layers  are  grown, 
a  small  amount  of  tellurium  is  added  to  the  solution.  The  melt  and  the  substrate  are  kept  apart 
during  the  first  part  of  the  run  by  a  small  ridge  in  the  bottom  of  the  boat.  Molten  stearic  acid 
is  then  poured  into  the  boat  to  cover  the  substrate  and  the  melt.  The  purpose  of  the  stearic 
acid  is  to  keep  the  surface  of  the  substrate  and  of  the  melt  clean  and  to  insure  good  uniform 
wetting  of  the  substrate  surface  by  the  melt.  The  boat  is  then  heated  by  the  graphite  heater 
strip.  At  about  300 °C  the  boat  is  tilted  so  that  the  melt  flows  over  the  substrate,  and  the  current 
in  the  heater  is  reduced  immediately  to  begin  cooling  the  boat.  When  the  melt  first  makes  contact 
with  the  substrate,  the  melt  is  not  quite  saturated  with  InSb  and  a  thin  layer  of  the  substrate 
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Fig.  1-2.  (a)  Apparatus  for  solution  regrowth  showing  mount  for  boron-nitride  boat, 
(b)  Inside  view  of  boron-nitride  boat. 


surface  goes  into  solution.  This  removes  any  surface  damage.  As  the  temperature  is  lowered, 
InSb  is  forced  out  of  the  solution  and  grows  in  a  uniform  layer  on  top  of  the  substrate.  At  about 
200°  to  250 °C  the  boat  is  tilted  back  and  the  remaining  solution  is  wiped  off  the  substrate.  The 
stearic  acid  is  then  poured  off  before  it  solidifies.  The  grown  surface  generally  has  a  clean 
mirror-like  appearance  and  the  growth  is  quite  uniform  across  the  substrate.  Figures  I-3(a-b) 
show  the  profile  of  a  grown  junction  delineated  by  etching  a  cleaved  (110)  surface.  All  our  solution 
regrowths  were  done  on  a  (1 00)  crystal  face  to  enable  cleaving  perpendicular  to  the  junction  plane. 

The  weights  of  the  various  ingredients  of  the  melt  used  in  growing  a  heavily  doped  n+  region 
on  a  p-type  substrate  area  of  12  x  12  mm  were  4.0  g  of  indium,  0.6  g  of  n-type  InSb,  and  4.0  mg 
of  tellurium.  The  weight  of  the  tellurium  dopant  in  this  case  was  chosen  to  produce  the  highest 
possible  doping  in  the  n  layer  (8  x  101  /cm  ).  '  The  ratio  of  the  InSb  weight  to  the  weight  of 
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indium  in  the  melt  was  chosen  from  published  solubility  data.  Figure  1-4  is  the  temperature 
variation  of  the  fractional  weight  of  InSb  required  to  produce  a  saturated  solution. 

A  typical  heating  cycle  is  shown  in  Fig.  1-5.  The  melt  is  poured  over  the  substrate  by  tilting 
the  boat  at  the  maximum  temperature  of  315°C  in  this  case.  At  about  225 °C  the  melt  is  wiped 
off.  Although  most  of  our  runs  were  similar  to  this  one,  the  times  and  temperatures  involved 
do  not  appear  to  be  critical. 


3.  Discussion 

The  main  innovation  of  our  solution-regrowth  technique  is  the  use  of  stearic  acid,  which 
facilitates  growing  at  the  low  temperatures  required  for  InSb.  Using  a  liquid  instead  of  a  gas 
atmosphere  also  greatly  simplifies  the  apparatus  and  the  process,  and  permits  visual  observa¬ 
tion  during  the  process.  Although  the  use  of  stearic  acid  is  limited  to  temperatures  not  higher 
than  350°C,  at  which  point  it  boils  away  too  rapidly,  other  liquids  such  as  boric  oxide  could  be 
used  for  crystals  which  require  higher  growth  temperatures  than  InSb. 
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Fig.  1-3.  (a)  Profile  of  120-p-thick  n  layer  grown  on  (100)  face  of  p-type  InSb  substrate.  Junction 

indicated  by  arrow  has  been  delineated  by  etching  cleaved  (110)  face,  (b)  Further  magnification  of 
profile  in  (a). 


Fig.  1-4.  Temperature  dependence  of  fractional 
weight  of  InSb  required  to  produce  saturated 
solution.  (After  Hall.^) 


Fig.  1-5.  Temperature  cycle  of  typical 
solution  regrowth. 
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One  application  of  this  method  was  the  fabrication  of  InSb  injection  lasers  in  which  the  light 
is  emitted  in  a  direction  normal  to  the  plane  of  the  junction.5  Such  a  structure  is  shown  in  Fig.  1-6. 
The  n+  layer  was  grown  on  high- resistivity  p-type  InSb  by  the  above  technique  and  the  p+  layer 
in  this  case  was  made  by  zinc  diffusion.  In  these  diodes  the  lifetimes  of  injected  carriers  are 
sufficiently  long  (10~7  to  lO’^sec)  and  the  mobilities  are  sufficiently  high  to  permit  population 
inversion  across  the  entire  region  between  the  contacts.  The  radiation  which  is  thus  generated 
in  the  active  region  has  to  pass  through  the  doped  p+  and  n+  layers;  hence,  these  have  to  be  of 
good  quality  crystal  which  is  transparent  to  the  radiation.  Solution  regrowth  proved  to  be  the 
only  satisfactory  technique  for  making  good  n+  layers. 

I.  Melngailis 
A.  R.  Calawa 


Fig.  1-6.  Large  volume  InSb  injection  laser  designed  for  coherent 
emission  normal  to  plane  of  junction. 

C.  INJECTION  MECHANISMS  IN  InSb  n+pp+  STRUCTURES 

We  have  studied  experimentally  the  injection  of  carriers  into  InSb  n*pp+  structures  similar 
to  those  used  as  bulk  lasers.^’ ^  In  particular,  we  have  determined  the  distribution  of  potential 
and  electric  field  along  the  p-region  at  various  current  levels  by  means  of  a  movable  potential 
probe  with  the  sample  at  77  °K.  These  measurements  indicate  that  at  currents  less  than  1  A/cm2 
the  dominating  mechanism  is  the  injection  of  electrons  into  the  p-region  from  the  n+  contact 
which  causes  a  saturation  of  electron  traps  and  an  increased  lifetime  in  the  p-type  bulk.  At 
higher  currents  an  appreciable  number  of  holes  is  injected  from  the  p+  region  in  order  to  neu¬ 
tralize  the  large  number  of  injected  electrons.  This  changeover  from  "single”  to  "double"  in¬ 
jection  is  observed  in  the  potential  distribution  as  well  as  in  the  patterns  of  light  emission  at 
the  surface  of  InSb  lasers. 

1.  Experimental  Technique 

Rectangular  bars  of  about  0.25  mm2  cross-sectional  area  and  lengths  varying  between  0.25 
and  3  mm  were  fabricated  from  p-type  InSb  with  a  resistivity  of  about  30  to  40  ohm-cm  at  77  °K. 
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The  p+  and  n+  contacts  were  alloyed  to  opposite  ends  of  each  bar;  in  some  cases  the  n+  region 
was  grown  by  a  solution-regrowth  technique.  A  bar  was  then  mounted  in  a  suitable  miniature 
package  and  connected  into  the  circuit  shown  in  Fig.  1-7.  Using  a  movable  tungsten  needle  probe 
to  measure  potential,  current- voltage  characteristics  could  be  obtained  from  any  number  of 
points  along  the  bar.  The  InSb  bar  and  tungsten  probe  were  immersed  in  liquid  nitrogen.  The 
probe  was  mounted  on  a  micromanipulator  and  was  viewed  through  a  calibrated  microscope. 
Sawtooth  pulses  0.1  to  2  msec  long  with  a  repetition  rate  of  about  100  pps  were  used  to  prevent 
heating. 

2.  Low  Current  Range 

At  low  currents  (0.5  to  1  A/cm  )  a  negative  resistance  region  is  observed  in  the  I-V  charac- 

8  9 

teristics  shown  in  Fig.  1-8.  This  negative  resistance  has  been  attributed  *  to  a  change  in  the 

lifetime  of  electrons  from  5  X  10”10  to  10~7  to  10-^sec  as  the  deep  traps  (0.05  eV  above  the 

10  + 
valence  band)  become  saturated  by  electrons  injected  from  the  n  contact.  The  traps  become 

saturated  first  near  the  n+  contact,  and  the  number  of  conduction  electrons  increases  sharply. 

As  the  current  is  increased,  the  electrons  penetrate  deeper  into  the  p- region  until  the  entire 
p-region  is  in  the  long-lifetime  high-conductivity  state.  The  potential  probe  measurements  of 
Figs.  1-8  and  1-9  illustrate  this  n propagation”  of  the  negative  resistance  along  the  bar  as  the 
current  is  increased.  We  observe  an  n+p  junction  potential  of  110  to  120  mV  (Fig.  1-9),  whereas 
the  p+p  contact  appears  ohmic  in  this  current  range.  This  is  consistent  with  the  electron  in¬ 
jection  model. 

Figures  I-lO(a-b)  show  the  potential  and  electric  field  distribution,  respectively,  in  the 
p-region  as  obtained  from  the  data  of  Fig.  1-8.  The  field  after  breakdown  (after  the  saturation 
of  traps)  is  small  (about  3  to  4  V/cm)  and  nearly  constant  with  distance  and  current.  The  field 
before  breakdown  is  nearly  constant  with  distance  but  increases  with  current  over  a  range  of 
about  5  to  70  V/cm.  The  data  of  Fig.  I-10(b)  show  that  just  before  breakdown  the  measured  field 
is  slightly  higher  than  that  expected  from  a  constant  resistivity.  A  possible  explanation  of  this 
is  a  variation  with  distance  of  the  quasi-Fermi  level  for  electrons  due  to  the  large  variation  in 
electron  density  near  the  boundary  of  the  plasma.  The  difference  in  quasi-Fermi  level  between 
two  points  with  electron  densities  n^  and  n^  is 


AV=  ~  log  — 
q  ee  n2 


Assuming  that  the  probe  measured  this  change  in  the  Fermi  level  for  electrons,  the  data  of 
Fig.  I- 10  yield  a  AV  (the  voltage  in  excess  of  the  expected  ohmic  voltage  drop)  of  about  80  to 
90  mV  over  a  distance  of  0.15  mm.  This  value  can  be  accounted  for  by  a  ratio  n^/n^  of  10^  to  10^. 


Fig.  1-7.  Circuit  used  to  measure  potential  as  a 
function  of  current  at  various  points  along  InSb 
bar.  Terminals  V  provide  the  horizontal  and 
terminals  I  the  vertical  input  to  oscilloscope. 
Source  consists  of  sawtooth  generator  and  am¬ 
plifier.  and  are  load  and  current  sam¬ 
pling  resistors,  respectively. 
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Fig.  1-8.  Potential  as  a  function  of  current  for  points  along  InSb  n  pp  device  146 
at  77° K.  Probe  settings  are  given  in  millimeters  measured  across  bar  from  n+p  junc¬ 
tion  to  p+p  junction.  The  two  traces  marked  n+  and  p+  were  obtained  with  probe 
just  inside  n+  and  p+  regions,  respectively. 


Fig.  1-9.  Potential  as  a  function  of  current  for  points  near  n+  junction  of  InSb  n+pp+ 
device  146  at  77°K.  Probe  settings  are  given  in  millimeters  measured  across  bar  from 
n+p  junction.  Trace  marked  n+  was  obtained  with  probe  just  inside  n+  region.  Slope 
of  this  n+  trace  is  due  to  resistance  between  n+  contact  and  mount,  and  it  is  to  be 
subtracted  from  slopes  of  remaining  traces. 
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Fig.  1-10.  (a)  Potential  and  (b)  electric  field  distribution  for  device  146 
at  77° K  showing  evidence  of  propagation  of  injected  plasma  across  bar 
with  increased  current.  Bar  length  is  approximately  1  mm. 
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This  large  a  ratio  maybe  possible  in  regions  where  the  lifetime  changes  rapidly  with  distance. 
There  is  also  an  added  potential  rise  in  regions  near  the  p+  contact  [as  seen  in  Fig.  I- 10(a)]  which 
cannot  be  explained  by  the  above  reasoning. 

3.  High  Current  Range 

2  -f 

At  currents  greater  than  10  A/cm  a  potential  begins  to  develop  across  the  p  p  junction. 

This  potential  increases  in  magnitude,  and  at  current  densities  of  100  to  200  A/cm2  it  becomes 
greater  than  the  potential  across  the  n+p  junction.  At  the  same  time,  a  low  field  region  (which 
increases  in  length  as  the  current  increases)  is  seen  near  the  p+  contact.  This  indicates  a  high 
carrier  density  near  the  p+  contact  resulting  from  the  injection  of  holes  from  the  p+  region. 
Figures  I-ll(a-b)  show  the  potential  and  field  distribution,  respectively,  for  an  n+pp+  InSb  device 
at  several  current  values  in  this  range.  The  low  field  regions  on  both  sides  of  the  bar  are  evi¬ 
dence  of  the  injection  of  carriers  from  both  contacts,  i.e.,  double  injection.  Throughout  the 
remainder  of  the  bar,  the  field  is  nearly  constant.  The  combined  length  of  the  two  low  field 
regions  in  the  0.5-mm  device  is  about  0.18  to  0.23  mm.  This  provides  an  estimate  for  the  opti¬ 
mum  length  of  an  InSb  bulk  laser,  since  the  largest  possible  active  volume  is  desired,  while 
avoiding  excessive  lattice  and  carrier  heating  which  can  occur  in  the  high  field  regions  of  the 
bar  in  Fig.  1-11 . 

Devices  were  also  fabricated  with  a  length  of  0.20  to  0.25  mm.  The  terminal  current- 
voltage  characteristic  for  one  of  these  devices  is  shown  in  Fig.  1-12.  The  potential  in  the  range 
between  1  and  500  mA  appears  mainly  across  the  two  junctions  and  the  current  increases  ex¬ 
ponentially  with  applied  voltage.  The  field  in  the  bulk  is  small  «3  V/cm)  and  remains  nearly 

constant  with  current.  This  implies  that  the  injected  carrier  density  increases  linearly  with 
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current  in  this  range.  However,  above  1  A  («  300  A/cm  )  the  field  in  the  bulk  begins  to  increase 
(15  to  25  V/cm  at  9  A). 

The  preceding  low  and  high  current  data  are  consistent  with  the  spontaneous  emission  pat¬ 
terns  observed  at  the  surface  of  InSb  bulk  lasers  at  11°K.  In  Fig.  1-13,  at  20  mA  the  brightest 
spot  and  hence  the  greatest  injected  carrier  density  is  near  the  n+  contact,  from  which  electrons 
are  injected.  At  50  mA  a  second  bright  spot  appears  near  the  p+  contact,  consistent  with  hole 
injection.  As  the  current  is  further  increased  to  250  mA,  the  light  intensity  and  hence  the  carrier 
density  near  the  p+  contact  increases  faster  than  the  density  near  the  n+  contact. 

R.  W.  Knepper 
I.  Melngailis 


D.  BULK  GaAs  MICROWAVE  AMPLIFIERS 

1 1 

Microwave  amplification  similar  to  that  recently  reported  by  Thim,  et  al.,  but  exhibiting 
some  significant  differences  in  frequency  and  threshold  conditions,  has  been  observed  in  samples 
of  n-GaAs  at  room  temperature. 

2 

All  the  samples  were  50-p-long  pieces  of  n-GaAs  with  a  mobility  of  about  5000  cm  /V-sec 
and  a  resistivity  of  approximately  lOOohm-cm.  Ohmic  contacts  were  made  to  opposite  sides  of 
the  sample  by  alloying  tin  at  400  °C  in  a  hydrogen  and  HC1  atmosphere.  The  samples  were 
mounted  in  a  varactor  diode  package  and  placed  at  the  end  of  a  50-ohm  coaxial  line.  The  reflec¬ 
tion  gain  was  then  measured  as  a  function  of  frequency  and  (pulsed)  bias  voltage.  The  gain 
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Fig.  1-11.  (a)  Potential  and  (b)  electric  field  distribution  for  device  153at77°K 
showing  carrier  injection  from  both  junctions.  Bar  length  is  approximately  0.5  mm. 
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Fig.  1-12.  Variation  of  current  with  total  voltage 
(including  junction  potentials)  for  device  162  at 
77° K.  Bar  length  is  0.  22  mm. 


20  mA  50  mA  250  mA 

Fig.  1-13.  Spontaneous  radiation  patterns  of  an  n+pp+  InSb  structure  with  400-p-long  bar-shaped 
p-region  obtained  at  11°K  with  three  different  values  of  current.  Diode  cross  section  was 
10“3  cm^.  Sensitivity  of  image  converter  was  reduced  as  current  was  increased. 
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Fig.  1-14.  Reflection  gain  as  a  function  of  frequency  for  an  applied  voltage  of  35  volts 
across  50-p-long  sample  of  n-GaAs,  whose  room-temperature  resistivity  was  100  ohm-cm. 


I-85-4614-1 1 


Fig.  1-15.  Current-voltage  characteristics  for  50-p-long 
samples  of  n-GaAs:  (a)  sample  whose  room-temperature 
resistivity  was  100  ohm-cm,  and  which  exhibited  gain  vs 
frequency  curve  shown  in  Fig.  1-14;  and  (b)  sample  whose 
room-temperature  resistivity  was  1  ohm-cm,  and  which 
did  not  exhibit  amplification. 
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remained  negative  (loss)  from  zero  applied  voltage  to  about  30  volts.  Above  30  volts,  gain  was 

observed  at  frequencies  near  3  GHz,  as  shown  in  Fig.  1-14  for  a  sample  biased  at  35  volts.  For 

biases  greater  than  36  volts  the  amplifier  output  became  noisy,  and  above  40  volts  the  sample 

oscillated  spontaneously  with  a  center  frequency  coincident  with  the  peak  of  the  gain  curve.  The 

saturation  output  power  of  this  sample  was  about  —  10  dBm,  consistent  with  the  dynamic  range 

reported  by  Thim,  et  al.,11  and  the  noise  figure  at  the  peak  of  the  gain  curve  was  about  17  dB 

for  10  dB  of  gain.  No  gain  was  observed  at  the  second  harmonic  of  the  gain  peak. 

The  frequency  at  which  the  gain  peak  occurs  agrees  fairly  well  with  the  frequency  at  which 
1 2 

Gunn  oscillations  occur  in  lower  resistivity  samples  of  the  same  length,  in  contrast  to  the  re¬ 
sults  of  Thim,  et  al.,  who  found  the  frequency  of  the  gain  peak  to  be  about  twice  the  Gunn  fre¬ 
quency.  Also,  the  average  applied  electric  field  necessary  for  gain  is  about  7000  V/cm—  twice 
the  value  reported  by  Thim,  et  al.  Both  these  differences  can  be  explained  if  the  electrical 
length  of  the  samples  used  by  Thim,  et  al.,  is  much  shorter  than  the  physical  length. 

In  Fig.  1-15,  the  current- voltage  curve  for  the  sample  of  Fig.  1-14  is  shown  and  compared 
with  that  for  a  sample  of  the  same  length  (50  p)  but  lower  resistivity  (1  ohm-cm)  which  exhibited 
oscillation  but  not  amplification.  We  see  that  the  amplifier  sample  exhibits  only  a  slight  tendency 
toward  current  saturation  with  increasing  voltage,  whereas  the  oscillator  sample  begins  to 
saturate  even  below  the  threshold  of  oscillation,  and  above  threshold  exhibits  a  range  of  current 
oscillation  relatively  independent  of  voltage. 

The  existence  of  a  positive  resistance  in  the  amplifier  samples  for  average  applied  electric 

fields  much  higher  than  the  field  at  which  Gunn  oscillation  normally  occurs  suggests  that  the 

usual  domain  formation  does  not  take  place  in  these  high- resistivity  samples.  This  agrees  with 

1 3 

recent  calculations  of  McCumber  and  Chynoweth  who  have  found  that,  for  sufficiently  small 

products  of  carrier  concentration  and  sample  length,  the  transferred  electron  effect  leads  to  a 

negative  conductance  at  frequencies  near  the  reciprocal  of  the  transit  time  even  though  the  DC 

14 

solution  displays  only  positive  resistance. 

A.  G.  Foyt 
T.  M.  Quist 
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II.  LASER  RESEARCH 


A.  NONCONVECTIVE  INSTABILITIES  IN  RAMAN  MEDIA 

An  investigation  of  the  dispersion  relation  for  infinite  homogeneous  Raman  media  repeals 
that  such  media  do  not  support  nonconvective  instabilities  near  the  Stokes  or  anti-Stokes  frequen¬ 
cies,  but  there  do  appear  to  be  such  unstable  modes  in  the  vicinity  of  cj  =  0,  k  =  0  and  cj  =  2cjt  , 
k  =  2k^,  where  oj^  and  k^  are  the  laser  frequency  and  wave  vector,  respectively.  The  time 
constant  for  the  growth  of  these  instabilities  is  on  the  order  of  milliseconds.  Because  the  as¬ 
sumed  dispersion  relation  describing  the  medium  is  not  reliable  for  these  values  of  oj  and  k, 
it  is  likely  that  these  instabilities  have  no  physical  significance. 

R.  S.  Sorbello 
H.  A.  Haus* 

B.  SELF-FOCUSING  OF  OPTICAL  BEAMS 


The  possibility  has  been  noted  of  the  self-trapping  of  optical  frequency  electromagnetic  beams 

due  to  a  nonlinear  index  of  refraction,  and  various  physical  mechanisms  which  could  give  rise  to 

1-3  3-5 

the  nonlinearity  have  been  discussed.  Trapped  modes  have  been  found.  Beam  trapping 
is  expected  to  have  important  effects  on  other  phenomena  occurring  in  intense  beams. 

Here  we  examine  the  possibility  of  an  increase  in  intensity  of  the  central  region  of  the  beam 
as  a  function  of  distance  in  the  direction  of  propagation  (self-focusing).  We  take  Eq.  (3)  of  Ref.  3 
as  our  starting  equation. 


v2?  ^ 

2  ,.2 
C  8t 


'  2  92(E2E) 

c2  at2  ' 


=  0 


(1) 


where  and  are  real,  and  e^E^  «  1.  We  assume  a  linearly  polarized  wave  of  frequency  c u 
and  propagating  along  the  z-axis,  and  obtain  the  approximate  equation 
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E'  =  0 


(2) 


where  ei,  =  3e^/4t  and  E'  is  the  slowly  varying  part  of  E  such  that  E  =  E'  exp  (ikz  —  ico t) .  Equa¬ 
tion  (2)  employs  a  familiar  approximation  of  optics  and  has  been  used  by  Vainshtein^  in  solving 
the  resonator  problem.  It  is  valid  when  the  distance  characteristic  of  the  change  of  E'  in  the  z 
direction  is  much  larger  than  a  wavelength.  The  equation  is  a  parabolic  type  and  can  be  solved 
with  an  open  boundary. 

From  Eq.  (2)  we  obtain  a  focusing  distance 

z  =  z-  =  i  a  JnA/n!,  /E1  (3) 

focus  2  v  O'  2  '  m 
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where  a  is  a  characteristic  transverse  radius  of  curvature  of  the  input  beam  intensity  and  E* 

m 

is  the  peak  field  value.  The  focusing  distance  is  defined  assuming  an  equiphase  input  beam.  It 
is  not  expected  that  the  solution  for  the  beam  center  will  be  strictly  correct  for  distances  within 
a  few  wavelengths  of  the  focusing  distance  because  of  the  neglect  of  the  second  derivative  inEq.  (2). 

To  verify  these  considerations,  Eq.  (2)  is  solved  numerically  for  a  cylindrically  symmetrical 
beam.  To  facilitate  this,  Eq.  (2)  is  rewritten  in  the  following  dimensionless  form 


.  9E*  92E*  1  9E * 

1  9  z*  +  9r*2  +  r*  97^ 


|E*  |2  E*  =  0 


(4) 


where  r*  =  r/a,  z*  =  z/2ka^,  and  E*  =  The  starred  notation  for  the  dimensionless 

quantities  is  similar  to  that  of  Ref.  3;  it  does  not  indicate  complex  conjugation.  In  this  dimension¬ 
less  form,  z?  =  (2  E*  )  *. 

focus  m  y 

We  have  solved  Eq.  (4)  by  a  finite  difference  or  mesh  technique  discussed  by  Harmutn  and 

modified  to  include  the  nonlinearity.  An  equiphase  Gaussian  intensity  profile  was  assumed  for 

the  input  beam.  An  example  of  the  numerical  results  obtained  on  an  SDS  930  computer  is  given 

in  Fig.  II- 1.  E*  was  chosen  large  enough  to  give  a  fairly  small  focusing  distance. 


Fig.  11-1.  Calculated  intensity  of  beam  center 

I*  vs  z*  (in  units  of  1/E*), 
m  m 
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The  curve  in  Fig.  II- 1  is  given  to  a  high  accuracy  by 

Im(z*)/Im(0)  =  z0(z02  “  Z*V1/2  (5) 

where  I  =  |E  |  ^.  We  find  that  z*  =  0.366/E*  (0),  while  from  Eq.  (3)  z*  =  0.353/E*  (0). 
m  1  m1  0  '  mv  M  focus  '  mv 

P.  L.  Kelley 


C.  HIGH  POWER  LASER  EFFECTS  IN  LIQUIDS 

Dielectric  breakdown  was  observed  in  water  by  focusing  a  ruby  laser  beam.  There  is  blue 
emission  from  the  focus,  less  prominent  than  in  air  breakdown  with  the  same  laser  output  and 
lens,  but  visible  through  a  blue  filter.  The  breakdown  has  a  threshold  at  about  1-MW  laser  power, 
where  the  blue  emission  changes  rapidly  with  slight  changes  of  laser  power.  At  lower  powers 
the  blue  emission  is  no  longer  detected,  either  visually  or  with  a  photomultiplier.  The  blue 
emission  has  a  half-duration  measured  as  20  nsec,  instrument  limited,  and  has  a  much  longer 
weak  afterglow.  Spectral  distribution  of  the  emission  was  measured  in  the  transverse  direction 

o  o 

from  3200  to  4500  A,  with  a  minimum  at  3600  A.  The  rise  toward  longer  wavelengths  corresponds 
to  the  region  of  the  Balmer  lines  of  hydrogen.  The  reason  for  the  relative  inconspicuousness  of 
the  blue  emission  may  be  that  the  Balmer  series  occurs  mainly  at  longer  wavelengths,  extending 

o 

to  6563  A.  Purification  of  the  water  shows  a  possible  increase  in  threshold  power,  but  does  not 
affect  the  spectral  distribution. 

Fluorescence  from  limonene  (C^H^)  was  induced  by  a  laser  beam,  presumably  by  multiple- 
photon  absorption,  and  observed  in  the  forward  direction.  In  the  wavelength  range  studied,  3200 

o  o 

to  4500  A,  fluorescence  occurs  as  a  band  from  4000  to  4500  A  and  beyond,  and  agrees  with  the 

normal  fluorescence  observed  on  a  spectrofluorometer.  The  horizontal  polarization  P^  is  about 

twice  the  vertical  polarization  P  ,  when  observed  with  incident  horizontally  polarized  light.  An 

8  ^ 

approximate  classical  theory,  describing  a  molecule  in  the  liquid  as  two  different  anisotropic 
oscillators,  one  polarized  by  the  incident  light  and  the  other  radiating  while  polarized  at  some 
other  orientation,  and  including  the  Brownian  rotation  of  the  excited  molecule  during  the  lifetime, 
gives  the  polarization 


P,  -  P 
h  v 

P,  +  P 
h  v 


r.  ,  1  x  kT  ,-1 

as  P=P0[1  +  (1-  3  PG)-^  r] 


(6) 


V  is  the  molecular  volume,  rj  the  viscosity,  r  the  lifetime  of  the  fluorescent  state,  and 

— 1/3^  p  ^l/2.  With  p  =  l/3  and  tj  =0.02  poises,  one  concludes  from  the  smallness  of  the  de- 
°  -9 

polarization  due  to  Brownian  motion  that  t  <  10  sec.  This  is  consistent  with  the  measurement, 
-8 

for  which  t  <  2  x  10~  sec,  instrument  limited. 

A  theoretical  examination  has  been  made  of  emission  of  sum  frequency  light  which  is  inco¬ 
herent  with  the  primary  light,  in  order  to  deduce  the  mechanisms  of  this  emission.  It  was  shown 
that  density  fluctuation  scattering  is  not  allowed  by  the  symmetry  of  the  molecular  polarizability 
tensors,  except  for  optically  active  molecules,  and  even  then  only  when  the  two  primary  frequen¬ 
cies  are  different.  In  other  cases,  therefore,  the  unshifted  sum  frequency  emission  must  arise 

from  orientation  fluctuations.  ~  T  .  , 

D.  L.  Weinberg 
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D.  HIGH  POWER  BRILLOUIN  AND  RAMAN  SCATTERING  IN  SOLIDS 

q 

After  completing  the  stimulated  Brillouin  scattering  experiment  in  quartz/  it  was  decided 
to  look  for  stimulated  Raman  scattering  in  quartz  in  order  to  study  possible  interference  between 
the  Raman  and  Brillouin  effects  at  high  incident  ruby  laser  powers.  Heretofore,  the  threshold 
for  stimulated  Raman  scattering  had  been  reached  in  only  a  very  few  solids  (calcite,  sulfur, 
diamond),  and  quartz  appeared  to  have  a  particularly  high  threshold. 

Figure  II- 2  shows  the  spectrum  of  quartz  photographed  in  the  forward  direction  (E  ||  c-axis) 
with  a  Bausch  and  Lomb  spectrograph  under  approximately  200-MW  ruby  laser  excitation  focused 
with  a  13-cm  lens  along  the  x-axis.  On  the  original  film,  four  stimulated  Stokes  and  two  anti- 
Stokes  lines  are  observable  with  the  following  wave -number  shifts. 


S  ^  =  459  cm”1 

AS ,  =  465  cm"1 
1 

-  921  cm"1 

AS^  =  930  cm-1 

=  1382  cm-1 

S.  =  1846cm”1 

-1 

All  measurements  are  to  ±1  cm  .  The  above  shifts  are  due  to  totally  symmetric  class-A  vibra¬ 
tions,  which  are  completely  polarized. 

The  difference  of  the  first  Stokes  shift  from  466  cm  1,  the  most  intense  line  in  the  quartz 

10 

spectrum  measured  by  classical  Raman  spectroscopy  methods,  seems  to  be  real.  The  spon- 

11  - i  _ 1 

taneous  line  is  at  least  6  cm  wide,  and  sometimes  the  weaker  452  cm  line  appears  on  the 

wings  of  the  466  cm  line.  Furthermore,  frequency  shifts  of  a  few  wave  numbers  of  some  of 

the  quartz  Raman  lines  have  been  noted  with  changes  in  relative  orientation  of  the  optic  axis  with 

respect  to  plane  of  scattering.1^  Finally,  first  overtones  observed  in  the  spontaneous  spectrum1  ^ 

_  \ 

fall  into  a  broad  band  consisting  of  harmonic  and  combination  lines  of  466  cm  and  its  satellites, 

-1  - 1  - 1 

while  second  overtones  are  not  observed  at  1398  cm  (3  x  466  cm  )  but  rather  at  1381  cm  ,  as 

in  the  present  experiment.  No  spontaneous  lines  have  ever  been  reported  in  the  vicinity  of 

- 1  - 1 

1864  cm  (4  X  466  cm  ),  whereas  spontaneous  AS^  cannot  be  observed  because  of  absence  of 
thermal  excitation. 


I- D0-3020-I I 


t  t  t  t  t  t  t 

as2  as^  laser  s<  s2  s3  s4 

Fig.  11-2.  Stimulated  Raman  spectrum  of  quartz  photographed  with  1.5-m 
Bausch  and  Lomb  spectrograph.  Ruby  laser  power  was  approximately  200  MW. 

An  interesting  property  of  quartz  is  that  the  dispersion  is  very  small  in  the  spectral  region 
around  the  ruby  laser  wavelength.  In  fact,  even  though  dn/dX  curves  slightly,  dn/dk  is  virtually 
flat.  This  means  that  phase  matching  is  determined  by  the  usually  small  secondary  term  in  the 
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phase -matching  condition,  and  AS^  is  observed  practically  in  the  forward  direction.  Since  the 
selection  rules  do  not  a  priori  forbid  scattering  from  one  kind  of  polarization  to  another,  it  might 
be  possible  to  observe  phase -matched  extraordinary  AS^  radiation  excited  from  incident  laser  radi¬ 
ation  polarized  as  an  ordinary  ray.  This  process  would  occur  at  a  relatively  large  scattering  angle. 

P.  E.  Tannenwald 
F.  H.  Perry 


E.  PHOTOELECTRON  STATISTICS  PRODUCED  BY  LASER 
OPERATING  BELOW  THRESHOLD  OF  OSCILLATION 

When  narrow -band  Gaussian  light  is  incident  upon  a  photomultiplier,  the  probability  distri¬ 
bution  of  the  photoelectrons  emitted  within  a  time  interval  short  compared  to  the  inverse  band- 

13- 17 

width  of  the  light  should  follow  the  Bose-Einstein  distribution.  Until  the  advent  of  the  laser, 

it  was  not  possible  to  make  counting  measurements  within  time  intervals  short  compared  to  the 
inverse  bandwidth  of  the  light.  A  laser  operating  slightly  below  threshold  may  emit  narrow-band 
light  with  a  Gaussian  amplitude  distribution  centered  about  one  frequency,  provided  one  axial 
mode  is  sufficiently  near  the  threshold  of  oscillation  so  that  its  gain  is  high  enough  compared  to 
the  other  modes  to  make  it  dominant.  By  "Gaussian  amplitude  distribution”  of  the  light  amplitude, 
we  mean  that  the  light  is  describable  in  Glauber's  P(a)  representation ^  in  which  each  mode 
has  a  P(a)  function  of  the  form 

P(a)  =  (l/7rn)  exp  (-|a  |  2/n)  (7) 

o 

For  the  experiments  to  be  reported  we  used  a  helium-neon  gas  laser  operating  at  6328  A  with  in¬ 
ternal  mirrors  of  hemispherical  geometry  and  reflectivities  of  greater  than  99  percent. 

By  monitoring  the  spectrum  of  the  photomultiplier  current  produced  by  the  light,  we  adjusted 
the  laser  to  operate  slightly  below  threshold  with  a  light  bandwidth  of  approximately  380 Hz.  In 
an  experimental  environment  the  vibrations  of  the  mirror  mounts  introduce  frequency  fluctuations 
which  cause  deviations  of  the  light  from  an  ideal  narrow-band  Gaussian  waveform.  Frequency 
fluctuations  are  not  detected  in  a  photomultiplier  homodyne  detector  and,  with  proper  care,  the 
influence  of  frequency  fluctuations  may  be  negligible,  as  demonstrated  experimentally  in  Ref.  18. 

It  is  legitimate,  therefore,  to  treat  the  waveform  emitted  from  the  laser  as  an  essentially  ideal 

19 

narrow-band  Gaussian  waveform.  The  power  output  of  the  laser  was  stabilized  against  long- 

_2 

term  drifts  by  a  feedback  network  with  a  time  constant  of  the  order  of  2  X  10  sec.  In  order  to 
reduce  background,  the  light  was  passed  through  a  polarizer  adjusted  for  maximum  transmission. 
Photoelectrons  emitted  in  the  photomultiplier  were  counted  in  a  high-speed  counter  with  a  double 

-9 

pulse  resolution  time  of  about  14  X  10  sec.  Three  sets  of  data  were  taken  for  counting  intervals 
-5  -4  -3 

of  10  ,10  ,  and  10  sec,  respectively.  The  number  of  photoelectrons  observed  in  successive 

time  intervals  was  recorded  on  tape  and  then  processed  by  a  computer.  Approximately  13,000 
samples  were  taken  for  each  of  the  three  counting  intervals  shown  in  Fig.  II- 3 .  If  the  laser  output 
indeed  consisted  of  a  pure  narrow-band  Gaussian  light  centered  at  a  single  frequency,  the  prob¬ 
ability  of  observing  exactly  n  photoelectrons  within  a  time  interval  short  compared  to  the  inverse 
bandwidth  would  follow  the  familiar  Bose-Einstein  Distribution  Law. 
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(a)  Counting  interval  T-=  10“5  sec; 
12,251  samples. 


(b)  Counting  interval  T  =  10“4  sec; 
13,619  samples. 


(c)  Counting  interval  T  =  10"3  sec; 
13,093  samples. 


Fig.  11-3. 


Probability  distribution  and 


number  of  observed  samples  vs  photoelectron  count. 
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P(n)  =  (1  +  n)"1  (1  4  l/nfn  (8) 

where  n  is  the  average  count.  On  the  semilog  plot  of  Fig.  II-3(a),  this  should  appear  as  a  straight 
line.  The  figure  shows  that  the  experimental  points  lie  on  or  about  a  straight  line  down  to  counts 
of  the  order  of  4.  At  lower  counts  the  experimentally  observed  points  deviate  from  the  theoret¬ 
ical  curve,  owing  to  admixture  of  other  modes  at  other  frequencies.  If  we  assume  that  these 
other  modes  are  broad-band,  and  therefore  effectively  produce  a  Poisson  probability  distribution 
of  their  own,  a  probability  distribution  slightly  different  from  the  Bose -Einstein  distribution  is 
calculated  to  be 


P(n)  =  (1  4  nm)_1  (1  4  l/nm)'n  exp  (nfa/nm) 


(9) 


where  n^  and  n^  are  the  average  photoelectron  counts  due  to  the  narrow-band  mode  and  the  broad¬ 
band  background,  respectively,  and  n^  =  n^  (1  +  l/n^).  Because  the  dominant  mode  is  linearly 
polarized,  the  background  can  be  estimated  from  the  degree  of  polarization  of  the  outcoming 
light.  Since  the  background  is  generally  unpolarized,  it  has  twice  the  intensity  of  the  light  po¬ 
larized  perpendicular  to  the  polarization  of  the  dominant  mode.  In  this  way  an  experimental  es¬ 
timate  of  the  background  can  be  made.  By  measuring  the  ratio  of  polarized  to  unpolarized  light, 

-5 

we  found  that  the  average  total  photoelectron  count  within  10  sec  was  15.8,  of  which  the  domi¬ 
nant  mode  and  the  background  counts  were  14  and  1.8,  respectively.  Introducing  these  values 
into  Eq.  (9)  results  in  the  theoretical  curve  shown  as  a  solid  line  in  Fig.  II-3(a).  Agreement  be¬ 
tween  this  theoretical  curve  and  the  experimental  results  is  evident.  The  experimental  points 
are  scattered  because  a  finite  number  of  samples  was  taken.  The  scattering  should  show  a 
mean-square  deviation  proportional  to  the  square  root  of  the  number  of  samples  taken.  The 
predicted  spread  is  indicated  by  vertical  bars  in  the  figures.  The  vertical  scale  given  for  S(n) 
on  the  right  side  of  Fig.II-3  indicates  the  number  of  samples  observed  for  any  given  photoelectron 
count  n.  For  high  values  of  n,  near  the  "tail"  end  of  the  figure,  the  number  of  samples  observed 
is  naturally  small  and  therefore  the  predicted  spread  indicated  by  the  vertical  bars  becomes  large. 
Also,  the  arrangement  of  all  experimental  points  along  horizontal  lines  which  correspond  to  the 
integer  values  of  S(n)  becomes  very  noticeable.  In  order  to  simplify  the  plotting  of  data  and  to 
reduce  the  scatter,  the  experimental  points  shown  in  Figs.  II -3 (b)  and  (c)  are  based  on  the  aver¬ 
ages  of  5  and  25  adjacent  data  points,  respectively.  Figure  II-3(b)  shows  the  results  obtained 

-4 

for  a  counting  interval  of  10  sec.  The  average  values  of  photoelectron  counts  used  in  Eq.  (9) 
to  compute  the  theoretical  curve  indicated  by  a  solid  line  in  Fig.II-3(b)  were  159,  140,  19  for 

_4 

the  total,  the  dominant  mode,  and  the  background,  respectively,  within  the  10  -sec  counting 

_  3 

interval.  Figure  II-3(c)  shows  the  corresponding  observations  for  a  counting  interval  T=  10  sec, 

-4 

which  is  longer  than  the  inverse  bandwidth  (in  radians)  of  the  light  (approximately  4.2  X  10  sec 
for  Fig.II-3).  In  this  range  an  asymptotic  expression  developed  by  R.  J.  Glauber  U  is  applicable, 

P(n)  =  (n/\l  (n)“3^2  exp  -  n/\/~n)2/2p]  (10) 
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where  \±  =  w/y,  y  is  the  linewidth,  w  is  the  average  rate  of  photoelectron  counts,  and  n  is  the 
average  photoelectron  count  within  the  counting  time  interval  T.  The  inverse  light  bandwidth  is 
determined  from  the  spectral  measurement  and  all  other  parameters  in  this  equation  are  ob¬ 
tained  from  the  photoelectron  counting  experiment.  Inserting  these  values  into  Eq.  (10),  we  ob¬ 
tain  the  theoretical  curve  shown  as  a  solid  line  in  Fig.  II-3(c).  The  shape  of  the  probability  dis¬ 
tribution  is  not  sensitive  to  small  admixtures  of  broad -band  modes  when  the  time  interval  be¬ 
comes  comparable  to  the  inverse  bandwidth.  Therefore,  the  equation  developed  by  Glauber  need 
not  be  corrected  for  the  background  modes.  Glauber's  equation  holds  for  counting  intervals  long 
compared  to  the  inverse  bandwidth.  In  the  present  example,  the  counting  interval  is  2.4  times 
longer  than  the  inverse  bandwidth,  yet  the  fit  of  Glauber's  formula  is  satisfactory. 

The  authors  believe  that  Fig.II-3(a)  shows  the  first  experimental  measurement  of  a  Bose- 
Einstein  probability  distribution  of  photoelectrons.  Also,  the  results  shown  in  Fig.  II-3(c)  con¬ 
stitute  the  first  experimental  test  of  Eq.  (10). 

A  determination  of  the  probability  distribution  of  photoelectron  counts  has  been  reported  by 

21 

Johnson,  McLean,  and  Pike,  using  an  incoherent  light  source  and,  therefore,  counting  inter¬ 
vals  much  longer  than  the  inverse  bandwidth.  Their  work  provided  the  initial  impetus  for  the 
measurements  reported  here.  ^  jrree(j 

H.  A.  Haus* 


*  Department  of  Electrical  Engineering  and  Research  Laboratory  of  Electronics,  M.I.T. 
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in.  MATERIALS  RESEARCH 


A.  Hg-Te  PHASE  DIAGRAM 

Differential  thermal  analysis  has  been  used  to  measure  liquidus  temperatures  for  the  por¬ 
tions  of  the  Hg-Te  system  from  3  to  25  and  from  50  to  100  atomic-percent  tellurium.  Liquidus 
points  between  30  and  60  atomic -percent  tellurium  were  obtained  previously1  from  data  on  the 
partial  pressures  of  Hg(g)  and  Te^fg)  in  equilibrium  with  Hg-Te  condensed  phases.  The  two  sets 
of  liquidus  points  have  been  combined  to  establish  the  phase  diagram  for  the  Hg-Te  system  shown 
in  Fig.  HI-1. 


Fig.  1 1 1-1 .  Phase  diagram  of  Hg-Te  system. 

The  melting  point  of  HgTe  measured  in  the  present  investigation  is  670  °C,  the  same  value 

2 

reported  by  Delves  and  Lewis.  The  temperature  and  composition  of  the  eutectic  between  HgTe 
and  tellurium  are  413.5  ±  0.5°C  and  84  atomic-percent  tellurium,  respectively,  in  comparison 
with  the  respective  values  of  411  ±  3°C  and  88  atomic-percent  tellurium  reported  by  Pellini  and 

-i 

Aureggi  in  1910.  For  samples  containing  80  or  more  atomic -percent  tellurium,  the  liquidus 
points  of  Pellini  and  Aureggi  are  in  good  agreement  with  the  present  values.  For  lower  per¬ 
centages  of  tellurium,  their  temperatures  are  systematically  lower  than  the  present  ones, 
probably  due  to  the  loss  of  mercury  from  their  samples  by  preferential  vaporization.  For  66 
atomic -percent  tellurium,  the  lowest  tellurium  content  they  studied,  their  liquidus  temperature 
is  548° C,  compared  with  the  present  value  of  583° C. 

Delves4  has  reported  a  region  of  two  immiscible  liquid  phases  in  the  vicinity  of  53  to  56 
atomic-percent  tellurium,  with  a  monotectic  temperature  of  664  ±  2°C.  The  present  data  do 
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not  confirm  the  existence  of  a  two-liquid-phase  region  in  this  location,  since  initial  liquidus  ar¬ 
rests  were  observed  at  661°  and  656°C,  respectively,  in  cooling  runs  on  samples  with  nominal 
compositions  of  54.2  and  55.0  atomic-percent  tellurium.  (The  actual  melt  compositions  at  the 
liquidus  temperatures  are  estimated  to  be  54.4  and  55.2  atomic -percent  tellurium,  respectively, 
because  of  preferential  vaporization  of  mercury.)  According  to  Delves'  results,  monotectic 
arrests  at  664  ±  2°C  should  have  been  observed  in  both  runs. 

A.J.  Strauss 
R.F.  Brebrick 


B.  PRESSURES  OF  Hg  AND  Se  IN  EQUILIBRIUM  WITH  HgSe(c) 

In  order  to  determine  the  pressures  of  mercury  and  selenium  in  the  vapor  over  Hg-saturated, 

Se-saturated,  and  congruently  subliming  HgSe(c)  between  450°  and  800° C,  the  optical  absorption 

°  5 

of  the  vapor  maintained  at  860°C  has  been  measured  between  2000  and  7000  A.  In  agreement 
with  mass  spectrographic  results,  no  evidence  for  the  species  HgSe(g)  was  found,  and  accordingly 
the  optical  density  was  assumed  to  be  the  sum  of  the  densities  due  to  mercury  and  selenium 
vapors.  The  partial  pressure  of  mercury,  and  the  sum  of  the  partial  pressures  of  all 

selenium  species,  p~,  could  therefore  be  obtained  from  the  measured  densities  together  with 

^1  6 
calibration  data  for  pure  mercury  and  pure  selenium. 

The  values  of  p^  and  p^  along  the  three-phase  lines  for  Hg-saturated  and  Se-saturated 

HgSe(c)  are  plotted  against  l/T  in  Fig.III-2.  Where  measured  values  (solid  lines)  could  not  be 


Fig.  MI-2.  Partial  pressure  of  Hg(g),  PHg'  an<^  sum  Par1‘*al 
pressures  of  all  selenium  species,  p^.,  along  three-phase  lines 
for  Hg-  and  Se-saturated  HgSe(c). 
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obtained,  the  plot  gives  calculated  values  (dashed  lines),  determined  for  one  component  from 
the  measured  values  for  the  other,  or  estimated  values  (dotted  lines).  Analysis  of  the  data 
shows  that  Hg-saturated  HgSe(c)  is  in  equilibrium  with  essentially  pure  Hg(f)  up  to  659°C  and 
contains  at  least  48.5  atomic -percent  selenium.  The  maximum  value  of  p^  over  HgSe(c)  is 
estimated  to  be  about  50  atm.  The  maximum  value  of  p^  is  0.62  atm,  the  pressure  over  Se- 
saturated  HgSe(c)  near  722°C. 

Values  of  p^  and  p^  in  equilibrium  with  congruently  subliming  HgSe(c)  were  measured  be¬ 
tween  about  450°  and  600°C.  The  data  are  well  represented  by  the  expressions  log  p^-g(atm)  = 

—  5.90(103)/T  +  6.22  and  log  p^fatm)  =  — 6.04(103)/T  +  5.82.  The  experimental  values  of  p^ 
agree  closely  with  those  calculated  from  the  measured  values  of  p^  by  applying  the  condition 
that  the  atom  fractions  of  mercury  and  selenium  in  the  congruently  subliming  vapor  are  essen- 

7 

tially  equal  and  using  the  published  equilibrium  constants  for  the  dissociation  of  Se4(g),  Se^(g), 
and  Seg(g)  into  Se2(g). 

Since  both  p^  and  pv  were  measured  for  congruently  subliming  and  Se-saturated  HgSe(c), 
these  data  could  be  used  to  calculate  the  standard  Gibbs  free  energy  of  formation  of  HgSe(c) 
according  to  the  expression  AG^.°[Hg(g)  +  l/2  Se2(g)  =  HgSe(c)]  =  RT  In  PpjgPz^2 •  The  values 
of  p^,  the  partial  pressure  of  Se2(g),  were  calculated  from  the  measured  values  of  p^,  by  using 
the  published7  equilibrium  constants  for  Se  vapor.  The  values  of  AG^°  obtained  for  congruently 
subliming  and  Se-saturated  HgSe(c)  are  the  same  within  the  limits  of  experimental  error.  Tak¬ 
ing  the  average  gives  AG°  =  —41.92  +  42.40(10  3)  T  kcal/mole  between  about  450°  and  600°C. 
These  are  the  first  values  which  have  been  reported  for  the  free  energy  of  formation  of  HgSe(c). 
Values  of  AH^°  =  — 10.8  kcal/mole  and  AS^°  =  —4.99  eu/mole  are  obtained  for  the  formation  of 
HgSe(c)  from  Hg(U  and  Se(c),  all  at  300°K,  by  extrapolating  the  free  energy  expression  and  using 
published  thermodynamic  data  for  HgSe(c),  Hg(g,U,  and  Se(g,  i>,c). 

R.F.  Brebrick 

C.  THERMODYNAMIC  STUDIES  OF  InTe  PRESSURE -TEMPERATURE  DIAGRAM 

g 

The  pressure-temperature  diagram  for  InTe  has  been  reported  previously.  The  pressure 
at  which  the  low-pressure  phase,  InTe  (I),  is  transformed  into  the  high-pressure  phase,  InTe  (II), 
decreases  with  increasing  temperature.  The  melting  point  of  InTe  (II)  increases  with  increasing 
pressure.  The  melting  curve  of  InTe  (I)  and  the  precise  location  of  the  InTe  (I)-InTe  (II) -InTe 
(liquid)  triple  point  could  not  be  determined.  It  was  found  that  at  sufficiently  low  temperatures 
InTe  (II)  can  be  retained  at  atmospheric  pressure  as  a  metastable  phase. 

As  an  extension  of  the  pressure-temperature  studies,  two  properties  of  InTe  have  been 
measured  at  atmospheric  pressure:  AV^.,  the  change  in  volume  of  InTe  (I)  on  melting,  and  AH^, 
the  heat  of  transformation  of  metastable  InTe  (II)  into  InTe  (I).  These  quantities  are  related  to 
the  slopes  of  the  melting  curve  of  InTe  (I)  and  of  the  phase  boundary  between  InTe  (I)  and  InTe  (II), 
respectively,  by  the  Clausius -Clapeyron  equation:  dT/dP  =  TAV/AH,  where  T  is  the  absolute 
temperature  of  the  phase  transition,  P  is  the  pressure,  AV  is  the  change  in  volume,  and  AH  is 
the  heat  of  transition. 

The  value  of  AVf  for  InTe  (I)  at  696°C  and  1  atm  was  determined  by  the  method  of  fed  and 
1  9 

unfed  castings,  as  modified  by  Ball.  The  densities  found  by  this  method  for  solid  and  liquid 
InTe  at  the  melting  point  are,  respectively,  5.845  ±  0.081  g/cm  (average  of  six  determinations) 
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Fig.  MI-3.  Recorder  traces  obtained  during  two  differential-scanning  calorimeter 
runs  in  which  high-pressure  form  of  InTe,  InTe  (II)/  was  transformed  into  low- 
pressure  form,  InTe  (I).  Abscissa  and  ordinate  give  temperature  in  degrees  Kelvin 
and  power  in  arbitrary  units,  respectively.  Dashed  lines  are  estimated  base  lines 
used  in  determining  areas  under  power  peaks,  (a)  Heating  rate:  20  deg/minute, 
(b)  Heating  rate:  40  deg/minute. 
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and  5.733  ±  0.033  g/cm  (average  of  seven  determinations).  The  corresponding  molar  volumes 
are  41.48  ±  0.44  and  42.28  ±  0.27  cm  /mole,  respectively,  and  the  value  of  AV  is  therefore  0.80 
±  0.72  cm^/mole.  Substituting  this  value  for  AV^.  and  the  measured10  value  of  4.29  kcal/g-atom 
for  AHj.  into  the  Clausius-Clapeyron  equation  gives  dT^/dP  =  2.l6deg/kbar  as  the  initial  slope  of 
the  melting  curve  for  InTe(I).  If  this  slope  does  not  change  with  pressure  up  to  the  triple  point, 
the  InTe(I)  melting  curve  intersects  the  InTe  (I)-InTe  (II)  phase  boundary  at  llkbar  and  718°C. 

g 

This  location  for  the  triple  point  is  close  to  the  one  estimated  previously. 

Values  of  AH^  for  the  InTe  (II)  -*■  InTe  (I)  transition  at  atmospheric  pressure  were  obtained 
by  both  metal-solution  and  differential-scanning  calorimetry.  The  first  method  was  used  to 
determine  AH^  at  0°C  by  comparing  the  temperature  changes  which  occurred  when  samples  of 
InTe  (I)  and  InTe  (II)  initially  at  0°C  were  separately  dissolved  in  molten  bismuth  initially  at 
3  50°C.  The  final  concentrations  of  indium  and  tellurium  in  the  bismuth  bath  were  kept  below 
2  atomic-percent  so  that  the  calculations  could  be  made  for  the  case  of  infinite  dilution.  The 

calorimeter  was  calibrated  by  adding  bismuth  initially  at  0°C  to  the  bismuth  bath.  Details  of 

1 1 

the  experimental  technique  and  calculation  procedure  are  given  elsewhere.  The  value  of  AH^ 
found  by  this  method  is  0.44  ±  0.01  kcal/g-atom. 

The  value  of  AH^  was  also  determined  with  a  differential-scanning  calorimeter  ( Perkin  - 
Elmer  Model  DSC-1).  In  this  apparatus  the  sample  and  an  inert  standard,  which  are  sealed  into 
separate  aluminum-foil  capsules,  are  heated  at  a  constant  rate  by  independent  resistance  heaters. 
The  differences  in  heater  power  required  to  keep  the  sample  and  standard  temperatures  equal 
throughout  the  run  are  recorded  as  a  function  of  temperature.  When  one  stable  phase  is  trans¬ 
formed  into  another,  a  peak  appears  on  the  recorder  chart  at  the  transformation  temperature, 
and  the  area  under  the  peak  is  proportional  to  AH^.  Although  there  is  no  thermodynamically 
defined  temperature  for  the  transformation  of  a  metastable  to  a  stable  phase,  AHt  for  such  a 
transition  can  also  be  measured  by  scanning  calorimetry  if  the  rate  of  transformation  is  strongly 
temperature  dependent.  In  that  event,  for  sufficiently  rapid  heating  rates  the  transformation 
occurs  predominantly  within  a  narrow  temperature  range  and  produces  a  well-defined  exothermic 
peak.  In  the  present  experiments,  two  such  peaks  were  observed  in  each  run,  as  seen  in  the 
typical  recorder  traces  shown  in  Figs.  III-3  (a-b).  These  traces  were  obtained  for  heating  rates 
of  20  and  40  deg/minute,  respectively. 

The  existence  of  two  exothermic  peaks  indicates  that  the  InTe  (II)  —  InTe  (I)  transformation 

12 

takes  place  in  two  steps,  and  therefore  confirms  the  finding  of  Sclar,  Carrison,  and  Schwartz 
that  at  atmospheric  pressure  an  intermediate  phase,  InTe  (II)’,  is  formed  during  the  transforma¬ 
tion.  In  eight  runs,  the  temperatures  at  which  the  evolution  of  heat  was  first  detected  were 

94  ±  6°C  for  the  first  peak  and  134  ±  6°C  for  the  second  peak.  These  values  are  consistent  with 

1 2 

the  results  which  Sclar,  et  al. ,  obtained  in  their  x-ray  and  metallographic  studies.  The  rela¬ 
tive  quantities  of  heat  liberated  during  the  two  steps  differed  considerably  from  run  to  run.  The 
percentage  of  the  total  evolved  during  the  first  step  ranged  from  14  to  39  percent.  However,  the 
total  heat  liberated  during  the  two  steps  was  quite  reproducible.  The  average  value  obtained  for 
AH^  is  0.42  ±  0.03  kcal/g-atom,  in  excellent  agreement  with  the  result  found  by  metal-solution 
calorimetry. 

The  value  of  AH  for  the  InTe  (II)  —  InTe  (I)  transformation  at  room  temperature  and  30  kbar 

t  g 

can  be  calculated  according  to  the  Clausius-Clapeyron  equation  from  the  measured  slope  of  the 
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Fig.  111-4.  Resistivity  at  77° K  (p 77)  vs  composition  for  the  high-pressure 
In Sb- In  system.  Composition  parameter  is  x  in  formula  In Sb  (11)^  (2  In)  . 


Fig .  II 1—5 .  Superconducting  transition  temperature  (Tc)  vs  composition  for  high- 
pressure  InSb-ln  system.  Where  there  are  two  bracketed  points,  two  discrete 
transitions  were  observed. 
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InTe  (I)-InTe  (II)  phase  boundary  and  the  value  of  AV  obtained  from  the  x-ray  densities  of  the 

t  13 

two  phases  at  atmospheric  pressure  and  their  compressibilities  between  1  atm  and  30  kbar. 

The  calculated  value  of  AH^  is  0.21  kcal/g-atom,  only  about  half  the  calorimetric  value  for  1  atm. 
The  difference  indicates  that  the  average  pressure  coefficient  of  heat  content  is  greater  for  InTe  (I) 
than  for  InTe  (II)  by  about  8  X  10  ^  kcal/(g-atom)-atm  between  1  atm  and  30  kbar. 

M.  D.  Banus  , 

P.  M.  Robinson’ 


D.  HIGH-PRESSURE  InSb-In  SYSTEM 

Additional  samples  of  nominal  composition  InSb(II)  (2  In)  have  been  prepared  in  order  to 

14  1_x  x 

extend  the  previous  investigation  of  the  high-pressure  InSb-In  system  over  the  entire  range 
from  pure  InSb(II),  the  high-pressure  phase  of  InSb,  to  pure  indium.  The  new  samples  were 

produced  at  37  kbar,  cooled  to  about  150°K  before  releasing  the  pressure,  and  stored  at  77°K, 

14 

as  described  previously,  in  order  to  retain  the  high-pressure  phase(s)  at  atmospheric  pressure. 

Resistivity  at  77°K  and  superconducting  transition  temperature  of  all  the  samples  are  plotted 

against  composition  in  Figs.  Ill -4  and  III-5,  respectively.  With  increasing  indium  content,  the 

resistivity  increases  to  a  maximum  at  about  x  =  0.3  and  then  decreases  smoothly  to  the  value 

for  pure  indium.  The  transition  temperature  initially  increases  sharply  from  the  value  of  2.1°K 

for  pure  InSb(II),  remains  essentially  constant  between  4.5°  and  5.0°K  from  x  =  0.1  to  x  =  0.7, 

and  then  decreases  to  the  value  for  pure  indium. 

14 

It  was  reported  previously  that  the  amount  of  heat  liberated  when  a  high-pressure  sample 
is  transformed  into  the  atmospheric  pressure  form  is  much  larger  for  x  =  0.2  than  for  x  =0.15 
or  below.  This  conclusion  was  based  on  time-temperature  data  obtained  in  drift  runs  during 
which  samples  were  heated  slowly  from  77° K  through  the  temperature  at  which  they  were  rapidly 
transformed  to  the  Stable  form.  In  additional  drift  runs,  it  has  been  found  that  the  amount  of 
heat  liberated  remains  relatively  high  between  x  =  0.2  and  x  =  0.5,  but  is  much  lower  for  x  =  0.6 
or  above.  We  conclude  from  this  result  that  there  is  a  phase  in  the  high-pressure  InSb-In  system 
which  differs  in  structure  from  both  InSb  (II)  and  indium.  This  phase  apparently  lies  between 
x  =  0.2  and  x  =  0.5,  but  its  exact  composition  cannot  be  determined  from  the  thermal  data. 

M.  D.  Banus 
Lynne  B.  Farrell 
A.  J.  Strauss 


E.  PRESSURE  DEPENDENCE  OF  B8  ~  B31  TRANSITION 
TEMPERATURE  IN  MnAs 

Under  atmospheric  pressure,  MnAs  is  transformed  at  44 °C  from  the  low-temperature  phase 
with  B8  structure  to  the  high-temperature  B31  phase.  Pressure  dependence  of  the  transition 
temperature  has  been  determined  by  measuring  the  change  of  resistance  with  hydrostatic  pres¬ 
sure  at  various  fixed  temperatures  between  25°  and  — 72°C.  The  experimental  technique  has 

15 

been  described  previously.  Results  of  the  measurements  are  shown  in  Fig.  III-6.  The  equa¬ 
tion  of  the  curve  drawn  through  the  data  points  is  P*  =  —5.6230  +  0.1073  T  —  0.000283  T^,  where 
P*  is  the  pressure  in  kilobars  and  T  is  the  transition  temperature  in  degrees  Kelvin.  The  form 
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Fig.  111-6.  Pressure  dependence  of  B8  ^  B31  transition 
temperature  in  MnAs. 


h2  x  io4  ( kG)2 


H  (kG) 


Fig.  111—7.  Magnetoresistance  (Ap/p0)  of  Ti203  as  a  function  of  magnetic 
field  (H)  and  of  H^.  Current  axis  is  tilted  45°  away  from  c-axis  toward 
a-axis.  Measurements  made  at  4. 2°K. 
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of  this  equation  suggests  that  there  is  a  critical  pressure  of  about  4.5kbar,  above  which  the 


B31  phase  is  stable  at  all  temperatures. 


J.  A.  Kafalas 
W.  A.  Newman 
J.  B.  Goodenough 


F.  ELECTRICAL  PROPERTIES  OF  Ti203 

Electrical  measurements  have  been  made  on  Ti^O^  single  crystals  containing  considerably 

less  nitrogen,  the  principal  impurity,  than  earlier  samples.  Reduction  in  nitrogen  content  has 

been  achieved  by  using  higher-purity  titanium  metal  in  preparing  the  oxide  and  by  using  hot 

titanium  as  a  getter  to  purify  the  argon  used  as  a  protective  atmosphere  in  synthesis  and  crystal 

16 

growth.  In  contrast  to  previous  results,  both  the  resistivity  (p)  and  Hall  coefficient  depend 
strongly  on  crystal  orientation,  and  p  now  increases  roughly  exponentially  with  decreasing 
temperature. 

Magnetoresistaiice,  Ap/p^,  has  been  measured  at  4.2°K  on  a  single  crystal  whose  c-axis  is 
oriented  at  45°  with  respect  to  the  current  direction.  As  shown  in  Fig.  III-7,  Ap/p^  is  propor¬ 
tional  to  the  square  of  the  magnetic  field  up  to  170  kG,  the  highest  field  employed.  The  value  of 
A p/pQ  at  170  kG  is  2.5,  compared  with  the  previous  maximum  value  of  0.12  for  Ti^O^  at  100  kG. 
The  value  of  2.5  is  too  high  to  be  explained  in  terms  of  a  one-band  model  for  the  transport 

Pr°PertieS*  J.M.  Honig 


G.  THEORY:  LOCALIZED  VS  COLLECTIVE  ELECTRONS 
IN  TRANSITION-METAL  COMPOUNDS 

There  are  two  limiting  theories  for  the  description  of  outer  electrons  in  solids:  a  localized- 
electron  (crystal-field)  theory  and  a  collective -electron  (band)  theory.  Whereas  f-electrons  are 
always  localized  and  outer  s-  and  p-electrons  are  almost  always  collective,  the  outer  d-electrons 
are  localized  in  some  crystals,  collective  in  others,  and  in  some,  localized  and  collective  d- 
electrons  are  simultaneously  present.  Spontaneous  crystallographic  distortions  of  transition- 
metal  compounds  from  a  higher  to  a  lower  symmetry  phase  can  provide  an  operational  criterion 
for  distinguishing  whether  the  outer  d-electrons  are  localized,  are  band  electrons  as  a  result  of 
cation-cation  interactions,  or  are  band  electrons  as  a  result  of  cation-anion  covalence. 

By  using  this  criterion,  it  is  possible  to  show  that  the  outer  d-electrons  in  CrS  are  localized, 
but  that  in  FeS  the  outer  d-electron  of  antiparallel  spin  is  collective,  belonging  to  the  cation  sub¬ 
lattice  as  a  whole  at  high  temperatures  and  to  molecular  clusters  at  low  temperatures.  The  TiP 
and  NbAs  structures  illustrate  changes  in  translational  symmetry  that  circumscribe  the  Fermi 
surface  of  cation-sublattice  band  electrons  with  a  Brillouin-zone  surface.  Distortions  from 
NiAs  to  MnP  structure  appear  to  be  induced  by  cation-anion  bonding,  cation-cation  interactions 
contributing  only  via  the  introduction  of  an  anisotropic  elastic  restoring  force. 

J.  B.  Goodenough 


H.  WET  CHEMICAL  ANALYSIS 

1.  Lanthanum  and  Yttrium  Cobalt  Oxides 

A  complexometric  method  accurate  to  three  parts  per  thousand  has  been  developed  for  the 
analysis  of  lanthanum  and  yttrium  cobalt  oxides.  The  samples  are  usually  dissolved  in  1:2  HNO^, 
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but  occasionally  it  is  necessary  to  fuse  the  sample  in  and  dissolve  the  fusion  in  10%  H^SO^. 

The  total  rare  earth  plus  cobalt  is  determined  on  an  aliquot  containing  about  10  mg  of  sample. 

The  aliquot  is  buffered  with  NH^NO^—  NH^OH  to  a  pH  between  8.3  and  8.6,  excess  EDTA  is  added, 

and  the  excess  is  titrated  with  standard  copper  solution  to  a  potentiometric  end  point  using  the 

17 

Reilley  gold-amalgam  electrode. 

The  rare  earth  alone  is  determined  on  an  aliquot  containing  about  5  mg  of  the  rare  earth. 

The  aliquot  is  buffered  as  for  tb^  above  determination,  the  solution  cooled  in  an  ice  bath,  and 
5%  KCN  solution  added  to  complex  the  cobalt.  After  two  minutes,  excess  EDTA  is  added.  Ex¬ 
cess  CN~  is  then  removed  by  adding  formaldehyde  dropwise  and  excess  EDTA  is  titrated  to  a 
potentiometric  end  point  using  the  gold-amalgam  electrode.  Automatic  titration  equipment  is 
used  for  both  titrations  of  excess  EDTA. 

2.  Iron-Nickel -Copper  Alloy  Thin  Films 

A  method  has  been  developed  for  the  analysis  of  iron-nickel-copper  alloy  thin  films.  Sam¬ 
ples  of  1  to  2  mg  are  dissolved  in  10%  HNO^.  Iron  is  determined  on  an  aliquot  containing  about 
0.2 mg  Fe  by  the  ortho-phenanthroline  spectrophotometric  method.  Nickel  interference,  which 
occurs  when  Ni/Fe  >  4,  is  eliminated  by  using  a  large  excess  of  ortho-phenanthroline  and  by 
making  the  photometric  measurement  immediately  after  the  ortho-phenanthroline  addition.  The 
spectrophotometric  measurements  are  made  in  a  1-cm  cell  at  495mjJL.  Beer's  law  is  obeyed  for 
solutions  containing  0.4  to  4.0  pg  Fe  per  milliliter.  The  method  is  accurate  to  about  three  parts 
per  thousand. 

Nickel  is  determined  by  the  usual  dimethylglyoxime  spectrophotometric  method.  Citric  acid 
is  added  to  eliminate  the  interferences.  Iodine  and  dimethylglyoxime  are  then  added  to  develop 
the  color,  which  is  measured  after  10  minutes  in  a  5-cm  cell  at  535  mjj..  Beer's  law  is  obeyed 
for  solutions  containing  0.4  to  2.0  pg  Ni  per  milliliter.  The  method  is  accurate  to  about  three 
parts  per  thousand. 

Copper  is  determined  on  an  aliquot  containing  0.4  to  2.0  jjl  Cu.  Citric  acid  and  EDTA  are 
added  to  eliminate  the  interferences,  the  pH  is  adjusted  to  8.5  to  9.0,  and  sodium  diethyl- 
dithiocarbamate  is  added  to  the  solution.  The  resulting  Cu  complex  is  extracted  in  CCl^,  and 
the  absorbance  of  the  CCl^  solution  is  measured  in  a  1-cm  cell  at  430  m\i.  As  suggested  by 
Sandell18  and  by  Snell,  Snell,  and  Snell,1  ^  the  aqueous  sodium  diethyl-dithiocarbamate  is  stored 
in  amber  glass  in  the  dark  until  use,  and  the  sample  solutions  are  protected  from  strong  light 
after  the  color  complex  is  formed.  This  system  does  not  obey  Beer's  law.  About  1 -percent 
accuracy  is  achieved.  E  B  0wens 

J.  C.  Cornwell 


34 


Section  III 


REFERENCES 


1.  R.  F.  Brebrick  and  A.  J.  Strauss,  J.  Phys.  Chem.  Solids  26,  989  (1965). 

2.  R.  T.  Delves  and  B.  Lewis,  J.  Phys.  Chem.  Solids  24,  549  (1963). 

3.  M.  Hansen  and  K.  Anderko,  Constitution  of  Binary  Alloys  (McGraw-Hill, 

New  York,  1958),  2nd  ed.,  p.  840. 

4.  R.T.  Delves,  Brit.  J.  Appl.  Phys.  16,  343  (1965). 

5.  R.  F.  Brebrick,  J.  Chem.  Phys.  43,  3846  (1965). 

6.  _ ,  J.  Chem.  Phys.  (in  press);  Solid  State  Research  Report, 

Lincoln  Laboratory,  M.  1.  T.  (1965:2),  p.36,  DDC  624611. 

7.  V.  V.  Illarionov  and  L.  M.  Lapina,  Doklady  Akad.  Nauk  S.  S.  S.  R.  114,  1021 
(1957). 

8.  M.  D.  Banus,  R.  E.  Hanneman,  K.  Gooen,  and  M.  Strongin,  Science  142, 

662  (1963). 

9.  D.  L.  Ball,  J.  Chem.  Eng.  Data  8,  61  (1963). 

10.  P.  M.  Robinson  and  M.  B.  Bever,  Trans.  Met.  Soc.  AIME  (in  press). 

11.  B.  W.  Howlett,  J.  S.D.  Leach,  L.  B.  Tickmar,  and  M.  B.  Bever,  Rev.  Sci. 

Instr.  33,  619  (1962). 

12.  C.  B.  Sclar,  L.  C.  Carrison,  and  C.  M.  Schwartz,  Science  147,  1569  (1965). 

13.  A.  J.  Darnell  and  W.  F.  Libby,  Phys.  Rev.  135,  A1453  (1964). 

14.  Solid  State  Research  Report,  Lincoln  Laboratory,  M.  I.  T.  (1965:2), 
p.  40,  DDC  624611. 

15.  Ibid.  (1965:1),  p.  50,  DDC  622446,  H-678. 

16.  Ibid.  (1965:2),  p.37,  DDC  624611. 

17.  C.N.  Reilley,  R.  W.  Schmidt,  and  D.  W.  Lamson,  Anal.  Chem.  30,  953  (1958). 

18.  E.  B.  Sandell,  Colorimetric  Determination  of  Traces  of  Metals  (Interscience, 
New  York,  1959),  p.437. 

19.  F.D.  Snell,  C.T.  Snell,  and  C.  A.  Snell,  Colorimetric  Methods  of  Analysis 
(Van  No  strand,  New  York,  1959),  0.46. 


35 


IV.  PHYSICS  OF  SOLIDS 


A.  ELECTRONIC  BAND  STRUCTURE 
1.  Reflectivity  of  ReOg 

In  the  last  Solid  State  Research  Report  (1965:2),  measurements  on  the  reflectivity  of  ReO^ 
were  reported  in  the  photon  energy  range  of  0.2  to  3.8  eV.  These  measurements  have  been 
extended  to  12  eV,  using  a  vacuum  ultraviolet  monochromator.  The  results  at  room  temperature 
are  shown  in  Fig.  IV- 1.  The  curve  shows  a  Drude-like  form  at  energies  less  than  2  eV  and  peaks 
at  4.6  and  9.6eV  due  to  interband  contributions.  By  approximating  the  reflectivity  for  low(<0.2eV) 
and  high  (>12.0  eV)  regions,  the  optical  constants,  n  and  k,  and  the  dielectric  constants,  e ^  and 

can  be  obtained  from  a  Kramers -Kronig  analysis  of  the  reflectivity  data.  In  theory,  the  die- 

c  f  f 

lectric  constant  e  =  +  ie^  can  be  separated  into  its  Drude  free-electron  parts  and  and  its 

bound-electron  parts  ej3  and  where 
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Fig.  IV- 1 .  Reflectivity  of  ReO^  at  room  temperature. 
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At  the  present  time,  various  techniques  are  being  used  to  separate  the  dielectric  constants 
and  obtained  from  the  Kramers-Kronig  analysis  into  free-  and  bound-electron  parts  in 
order  to  evaluate  such  parameters  as  m*  and  y  which  may  be  compared  with  DC  conductivity 
measurements.  An  attempt  will  also  be  made  to  interpret,  from  energy  band  theory  considera¬ 
tions,  the  origin  of  the  interband  transitions  above  2eV. 

W.  J.  Scouler  A.  Ferretti 
J.  Feinleib  B.  Feldman 


/m*  and  y  is  the  relaxation  frequency  which  is  not  necessarily  a  constant. 


2.  Band  Structure  of  ReOg 

1 

The  simplified  LCAO  (linear  combination  of  atomic  orbitals)  or  tight  binding  scheme  has 
been  used  to  investigate  the  band  structure  of  ReO^.  The  study  illustrates  the  connection  between 
ligand  field  theory  and  energy  band  theory,  and  indicates  the  range  of  validity  of  the  former  model 
for  this  compound.  Using  a  linear  combination  of  rhenium  6s-,  6p-,  and  5d-orbitals  and  of  oxygen 
2s-  and  2p-orbitals,  the  21  X  21  secular  determinant  was  solved  analytically  under  certain  sim¬ 
plifying  assumptions.  In  conformity  with  the  scheme  proposed  by  Goodenough,2  the  energy  levels 
were  found  to  separate  into  a  and  7r  bonding  and  antibonding  bands  and  into  a  set  of  nonbonding 
7r  levels  made  up  of  oxygen  p-orbitals.  Analytic  expressions  have  been  found  for  the  structure 
of  the  antibonding  7r  band  in  which  the  Fermi  level  of  ReO^  is  situated.  In  order  to  obtain  a 
quantitative  description  of  the  energy  bands,  estimates  of  the  various  parameters  which  occur 
in  the  simplified  LCAO  scheme  must  be  made.  This  is  being  accomplished  by  using  the  augmented 
plane-wave  (APW)  method  of  Slater  to  fit  these  bands  to  electronic  energy  values  calculated  at 
symmetry  points  T  and  X,  and  along  A,  the  line  joining  these  points.  This  permits  detailed 
predictions  concerning  the  magnetic  and  transport  properties  of  ReO^. 

J.  M.  Honig 
J.  O.  Dimmock 


3.  Magnetoabsorption  of  Direct  Transition  in  Gallium  Selenide^ 

Investigation  of  the  band  structure  of  GaSe,  by  studying  interband  magneto-optical  absorp- 

4 

tion,  is  being  continued.  Much  more  highly  resolved  data  than  that  reported  previously  have 
been  obtained  in  the  Faraday  geometry.  Preliminary  results  have  also  been  obtained  in  the 


tThis  experiment  was  conducted  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.T. 
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Voigt  configuration  (magnetic  field  in  the  plane  of  the  layers  and  the  propagation  vector  of  the 
radiation  along  the  c-axis);  the  latter  observations  were  made  at  1.5  °K  and  with  magnetic  fields 
of  zero  through  90  kG. 


Fig.  IV-2.  Optical  absorption  near  direct  edge 
in  GaSe  at  1 ,5°K  with  propagation  vector  T  |  | 
c-axis. 


Samples  have  now  been  studied  for  which  the  zero-field  data  at  1.5  °K  show  as  many  as  three 
clearly  defined  exciton  peaks  in  addition  to  the  well-known  broad  absorption  peak  above  the  con¬ 
duction  band  edge.  A  plot  of  the  zero-field  absorption  constant  as  a  function  of  energy  for  one  of 
the  better  samples  is  given  in  Fig.  IV-2.  The  binding  energy  of  the  exciton  as  determined  from 
this  series  of  lines  is  (21  ±  1)  meV.  Assuming  that  GaSe  is  essentially  a  three-dimensional 
material  and  using  a  dielectric  constant  of  10.6,  which  is  obtained  from  reflectivity  data,  one 
calculates  a  reduced  effective  mass  m*  =  0.17  m  from  the  observed  exciton  binding  energy  of 
21  meV. 

A  plot  summarizing  some  of  the  more  recent  magnetoabsorption  data  taken  in  the  Faraday 
geometry  is  shown  in  Fig.IV-3.  The  data  shown  are  for  left  circularly  polarized  light;  data 
taken  for  right  circularly  polarized  light  give  a  similar  spectrum  which  is  very  slightly  displaced 
to  higher  energies.  The  triangles  represent  data  obtained  by  sweeping  the  wavelength  at  a  fixed 
value  of  magnetic  field;  the  crosses  represent  data  obtained  by  sweeping  the  magnetic  field  at  a 
fixed  value  of  wavelength.  Magnetic  quantum  numbers  of  the  particular  Landau  sub-bands  cor¬ 
responding  to  the  observed  transitions  have  been  identified  and  are  indicated  in  the  figure.  The 
transitions  themselves  are  all  actually  to  exciton  levels  lying  beneath  these  sub-bands.  From 
the  higher  quantum  number  lines,  a  reduced  effective  mass  of  m*  =  0.14m  is  obtained.  From 
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Fig.  IV— 3,  Summary  of  magnetoabsorption  results  in 
GaSe  at  1 ,5°K  with  magnetic  field  H  1 1  r||  c-axis 
and  left  circular  polarization.  Spectrum  for  right 
circular  polarization  is  slightly  displaced  and  has 
been  omitted  for  clarity. 


Fig.  IV-4.  Summary  of  magnetoabsorption  results  in  GaSe 
at  1 .5°K  with  Hi  c-axis  and  Tllc-axis. 
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the  convergence  of  the  higher  n-number  "quasi-Landau  level"  lines,  a  zero-field  edge  energy 
of  2.130  eV  and  an  exciton  binding  energy  of  -20  meV  are  obtained;  this  effective  mass  and  bind¬ 
ing  energy  is  in  agreement  with  the  values  obtained  from  the  zero-field  exciton  spectrum. 

The  magnetoabsorption  spectra  in  the  Voigt  configuration  at  1.5  °K  show  oscillatory  effects 
which  are  appreciably  smaller  in  amplitude  and  more  complex  than  those  obtained  in  the  Faraday 
geometry.  A  plot  of  the  energies  corresponding  to  the  transmission  minima  as  a  function  of 
magnetic  field,  as  shown  in  Fig.  IV-4,  gives  a  collection  of  lines  of  nonuniform  spacing  and  ap¬ 
pears  to  indicate  the  presence  of  more  than  one  series  of  transitions.  In  the  present  preliminary 
experiments  in  the  Voigt  geometry,  we  have  not  as  yet  carried  out  studies  for  different  orienta¬ 
tions  of  the  magnetic  field  in  the  c-plane. 

There  are  several  models  for  the  band  structure  of  GaSe  which  are  consistent  with  the 
present  experimental  data.  The  three  most  likely  possibilities  are  (a)  simple  valence  and  con¬ 
duction  bands,  both  located  at  the  T-point  in  the  zone,  (b)  cigar-shaped  ellipsoidal  valence  and 
conduction  bands,  located  at  the  three  equivalent  vertices  of  the  hexagonal  Brillouin  zone  (Z- 
point),  or  (c)  a  cigar-shaped  ellipsoid  of  revolution,  located  at  T,  for  either  the  valence  or  con¬ 
duction  band  and  a  more  complex  form  for  the  remaining  band. 

J.  Halpern 


4.  Cyclotron  Resonance  in  Gallium  Selenide 

Measurements  of  the  microwave  conductivity  of  GaSe  as  a  function  of  magnetic  field  are 
being  carried  out  with  the  hope  of  observing  cyclotron  resonance.  This  solid  is  of  particular 
theoretical  interest  because  of  its  layered  structure.  The  large  oscillatory  magnetoabsorption 
effects  obtained  by  Halpern  (Sec.  IV-A-3)  on  similar  material  indicate  that  the  carrier  relaxation 
times  may  be  sufficiently  long  to  resolve  cyclotron  resonance. 

In  the  experiment,  the  sample  is  cemented  to  the  bottom  surface  of  a  TMQ13  cylindrical 
cavity  resonant  at  70  GHz  and  cooled  to  4.2°K.  The  static  magnetic  field  is  applied  in  the  plane 
perpendicular  to  the  c-axis  and  is  always  perpendicular  to  the  RF  electric  field.  Carriers  are 
excited  by  illuminating  the  sample  with  tungsten  light  chopped  at  260  Hz  and  the  signal  reflected 
from  the  cavity  is  phase-sensitive  detected  at  the  same  frequency. 

Out  of  six  samples  examined  so  far,  one  has  shown  a  well-resolved  absorption  peak  at  a 
field  of  500  G,  corresponding  to  an  effective  mass,  m*  =  0.02  m,  and  with  a  line  shape  that  in¬ 
dicates  a  value  of  approximately  1.5  for  a ;r.  No  significant  shift  in  the  position  of  the  peak  is 
observed  as  the  magnetic  field  is  rotated  in  the  plane  perpendicular  to  the  c-axis.  Only  one 
peak  is  observed  in  fields  up  to  25  kG  and  it  is  not  known  whether  this  corresponds  to  holes  or 
electrons.  Resonance  with  the  magnetic  field  in  a  plane  containing  the  c-axis  has  not  yet  been 
investigated. 

Because  the  above  result  has  been  observed  in  one  sample  only,  it  must  be  regarded  as 
preliminary  and  not  conclusive.  Efforts  are  continuing  to  observe  effects  in  other  samples. 

S.  Bermon 
W.  C.  Kernan 
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5.  Magnetoplasma  Cyclotron  Resonance  in  p-Type  Lead  Selenide 

Measurements  of  the  microwave  surface  impedance  at  70  GHz  as  a  function  of  magnetic  field 

18  3 

of  single  crystal  p-type  PbSe  with  a  carrier  concentration  of  10  holes/cm  are  continuing.  The 
conditions  of  the  experiment  are  such  that  one  does  not  expect  to  observe  the  usual  cyclotron 
resonance,  as  for  example  in  germanium,  but  resonances  strongly  affected  by  the  high  density 
plasma  which  exists  in  the  solid. 

In  the  present  experiment  the  static  magnetic  field  is  in  the  plane  of  the  sample  surface  which 
forms  the  central  section  of  the  bottom  wall  of  a  cylindrical  cavity  operating  in  the  TE113  mode. 
The  sample  surface  itself  must  be  of  an  extremely  high  quality,  free  from  strain  or  damage,  to 
satisfy  the  requirements  on  the  relaxation  time,  cur  >  1.  A  mica  strip  placed  over  the  bottom 
surface,  along  a  diameter  of  the  cavity,  removes  the  degeneracy  of  the  cavity  mode  and  allows 
the  direction  of  the  RF  electric  field  to  be  fixed.  Experiments  are  carried  out  at  4.2°  and  1.2  °K. 

Absorption  derivative  curves  in  the  range  below  3  kG  exhibit  marked  structure  with  very 
good  repeatability.  Using  the  classical  skin  effect  theory  and  an  infinite  relaxation  time,  a  pre¬ 
liminary  theoretical  analysis  has  been  carried  out  for  the  data  obtained  with  the  static  magnetic 
field  aligned  along  the  [100]  and  [110]  axes,  and  with  the  ordinary  and  extraordinary  modes  of 
propagation  in  the  crystal  excited  independently.  Points  of  sharp  descent  are  identified  with 
dielectric  anomalies.  The  behavior  of  the  structure  is  consistent  with  that  predicted  by  using 
the  values  of  effective  mass  and  anisotropy  obtained  from  previous  Shubnikov-de  Haas  measure- 

5 

ments  on  this  sample.  Azbel-Kaner-type  resonance  as  an  explanation  for  the  main  structure  is 
ruled  out.  A  more  complete  theoretical  treatment  which  includes  the  effect  of  a  finite  relaxation 
time  requires  machine  calculation.  Such  a  program  is  currently  under  way. 

Additional  weak  structure,  which  appears  to  resemble  Shubnikov-de  Haas  oscillations,  has 
also  been  observed  at  high  magnetic  fields.  These  effects  are  being  examined  further. 

S.  Bermon 
W.  C.  Kernan 


6.  Fourier  Expansion  for  Energy  Bands  in  a  Periodic  Solid 

Electronic  energy  bands  in  solids  are  periodic  in  wave-vector  space  and  hence  are  amenable 
to  a  Fourier  expansion.  The  symmetry  of  the  lattice  and  the  degeneracies  of  the  bands  at  high 
symmetry  points  in  the  Brillouin  zone  are  included  by  using  a  generalization  of  the  technique 
developed  by  Luttinger^  for  obtaining  a  Taylor's  series  expansion  of  the  energy  bands.  The 
effective  mass  Hamiltonian  thus  derived  is  a  matrix  with  dimension  equal  to  the  number  of  bands 
being  considered.  Experimentally  observed  quantities  such  as  cyclotron  effective  masses,  Fermi 
surface  cross  sections,  de  Haas-van  Alphen  periods,  etc.,  are  readily  computed  and  thus  can  be 
used  experimentally  to  determine  the  Fourier  coefficients  of  the  energy  bands. 

This  technique  is  being  applied  to  obtain,  phenomenologically,  the  energy  bands  in  the  semi¬ 
metal  antimony,  where  recent  magnetoreflection  experiments^  have  given  information  about 
energy  band  gaps  and  effective  masses  at  points  L  and  T  in  the  Brillouin  zone.  In  order  to 
establish  the  feasibility  of  such  expansion  in  a  metal,  an  application  to  aluminum  is  also  being 

undertaken.  ^  ^  , , 

G.  F.  Dresselhaus 

Mildred  S.  Dresselhaus 
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B.  HYPERSONIC  WAVES  IN  SOLIDS 

1.  Generation  of  Hypersonic  Waves  by  Using  Thin-Film  CdS  Transducers 

Many  of  the  difficulties  and  limitations  arising  from  the  excitation  of  hypersonic  waves  by 
surface  effects,  harmonic  quartz  transducers,  or  magnetostrictive  films  may  be  avoided  by 
using  vapor -deposited  insulating  films  of  piezoelectric  materials.  Consequently,  an  electron- 
beam  technique,  to  be  described  in  a  later  report,  has  been  developed  for  the  vapor  deposition 
of  insulating,  piezoelectric  thin-film  CdS  transducers  onto  various  substrates. 

The  ideal  CdS  film  would  be  a  single  hexagonal  crystal  with  its  c-axis  perpendicular  to  the 
film  surface.  Then,  by  orienting  the  driving  electric  field  with  respect  to  the  film  surface,  it 
would  be  possible  to  excite  either  longitudinal  or  transverse  hypersonic  waves;  these  waves 
could  then  be  used  as  a  probe  to  study  certain  properties  of  the  substrate  material. 

To  date,  acoustically  active  CdS  films  have  been  deposited  on  single  crystal  specimens  of 
aluminum,  permalloy,  indium  antimonide,  gallium  arsenide,  quartz,  sapphire,  and  ruby.  Pulse- 
echo  techniques  were  used  to  determine  whether  or  not  a  given  film  was  acoustically  active. 

The  electromagnetic  frequency  range  of  excitation  of  the  longitudinal  mode  in  active  films  varied 
from  30.2  MHz  for  CdS  on  an  aluminum  substrate,  to  70  GHz  for  CdS  on  a  z-cut  quartz  rod 
(Sec.  IV-B-2).  The  corresponding  thickness  range  was  from  70  to  less  than  0.1  jjl . 

Whenever  films  thicker  than  one -half  a  wavelength  have  been  deposited  on  quartz,  transduc¬ 
tion  efficiency  has  improved  as  compared  with  surface  generation  in  quartz  alone.  Figure  IV-5(a) 
indicates  the  echo  pattern  obtained  in  the  transverse  mode  at  1  GHz  with  a  0.9-p  film  deposited 
on  an  x-cut  quartz  substrate.  Figure  IV-5(b)  shows  a  similar  pattern  obtained  at  1  GHz  for  the 
same  substrate  but  without  a  film. 

The  patterns  shown  in  Fig.IV-6(a)  were  obtained  for  propagation  in  the  longitudinal  mode 
at  1GHz  through  z-cut  quartz  without  a  deposited  film  and  those  in  Fig.  IV-6(b)  with  a  2.2 -p  film 
deposited  on  the  same  piece  of  quartz  under  identical  conditions. 


Fig.  IV-5.  Pulse-echo  pattern  of  1-GHz  transverse 
waves  in  x-cut  quartz  at  300°K/  (a)  generated  and 
detected  by  CdS  transducer  0.9p  thick  and  (b)  gen¬ 
erated  and  detected  by  surface  generation  in  quartz 
alone. 
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Fig.  IV-6.  (a)  Expected  result  of  attempt  to 

excite  1-GHz  longitudinal  waves  in  z-cut  quartz 
at  room  temperature  by  surface  excitation  of 
quartz,  (b)  Pulse-echo  pattern  of  longitudinal 
waves  in  this  same  specimen  of  z-cut  quartz  under 
identical  conditions  but  with  vapor-deposited  CdS 
transducer  2.2  p  thick  generating  and  detecting 
hypersonic  waves. 


Fig.  IV— 7.  Pulse-echo  pattern  of  longitudinal 
waves  in  polycrystalline  permalloy,  (a)  at  41  MHz 
and  300°K,  generated  and  detected  by  vapor- 
deposited  CdS  transducer  of  thickness  5.4p  and 
(b)  at  40.  8  MHz  and  300° K,  with  the  waves  gen¬ 
erated  and  detected  by  resonant  quartz  plate 
bonded  to  sample  and  operating  in  its  third 
harmonic. 


Fig.  IV-8.  "Variation  of  asymptotic  Curie  temperature 
with  Cr3+-Cr3+  separation. 
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Figure  IV-7  shows  the  echo  patterns  for  the  longitudinal  mode  propagation  in  a  permalloy 
sample,  comparing  a  directly  deposited  CdS  film,  with  a  fundamental  frequency  at  41  MHz 
[Fig.  IV-7(a)],  to  a  resonant  quartz  plate  transducer  bonded  to  the  same  piece  of  permalloy  and 
operating  in  its  third  harmonic  at  40.8  MHz  [Fig.  IV-7(b)]. 

Immediate  future  plans  call  for  measuring  the  transduction  efficiency  of  CdS  films  as  a 
function  of  frequency  and  of  mode  of  generation.  A  parameter  of  interest  will  be  the  substrate 
material.  Concurrently,  an  x-ray  analysis  will  be  attempted  in  order  to  determine  the  stoichi¬ 
ometry  of  the  films  and  its  correlation  with  the  piezoelectric  efficiencies  for  the  different  modes 
of  generation.  R  Weber 

J.  W.  Burke 


2.  Microwave  Phonon  Generation 

o 

Longitudinal  phonon  echoes  have  been  observed  for  the  first  time  at  70  GHz  using  1600 -A 

thick  CdS  film  transducers,  which  were  deposited  on  quartz  by  R.  Weber  and  J.  W.  Burke 

5 

(Sec.  IV-B-1).  The  measured  velocity  in  quartz  of  6.25  x  10  cm/sec  agrees  within  experimental 

error  with  the  value  reported  in  the  literature  for  the  velocity  of  the  longitudinal  mode  in  the  z- 

direction.  Four  echoes  were  observed,  with  the  first  echo  6dB  above  noise  level.  Although  the 

overall  efficiency  of  the  transduction  process  for  this  particular  quartz  sample  and  CdS  film  is 

less  than  that  usually  observed  for  x-cut  quartz,  no  general  statement  can  be  made  about  the 

performance  of  CdS  films  without  testing  at  least  several  samples.  Successful  development  of 

this  technique  will  allow  an  investigation  of  the  propagation  of  sound  in  z-cut  quartz,  c-axis 

sapphire,  and  other  nonpiezoelectric  materials.  T 

J .  Id.  1  naxter 

C.  D.  Parker 


C.  MAGNETISM 

1.  Magnetic  Properties  of  Chromium  Spinels  with  Nonmagnetic  A-Sites 

Study  of  the  magnetic  properties  of  the  spinels  ACr^X^,  where  A  is  nonmagnetic  zinc  or 
cadmium  and  X  is  oxygen,  sulfur,  or  selenium,  first  reported  in  the  previous  Solid  State  Re- 

g 

search  Report,  has  been  continued.  We  have  found  that  CdCr2S4,  like  CdCr2Se4,  is  ferromag¬ 
netic  with  a  moment  of  5.55  at  4.2 °K  in  10,000  oe,  a  Curie  temperature  T^  =  86 °K,  a  Curie- 
Weiss  constant  C^  =  3.75,  and  an  asymptotic  Curie  point  T^  =  156 °K.  We  had  found  that 
CdCr204  becomes  antiferromagnetic  below  T^  =  9°K.^  Measurements  of  the  paramagnetic  be¬ 
havior  show  that  Cjyj  =  3.  54  and  T^  =  —  88  °K  in  this  material. 

The  effect  of  cell  size  on  T^  in  the  entire  series  of  compounds  is  shown  in  Fig.  IV-8.  It 
can  be  seen  that  there  is  a  strong  dependence  of  T^  on  lattice  size  for  either  A-site  ion,  but  that 
the  quantitative  variation  is  quite  different  depending  on  whether  the  A-site  ion  is  zinc  or  cadmium. 
Since  T^  is  proportional  to  the  sum  of  all  the  magnetic  exchange  interactions,  this  further  indi¬ 
cates  that  the  diamagnetic  A-site  cation  significantly  affects  these  interactions.  Specifically,  the 
data  shown  in  Fig.  IV-8  indicate  that  for  a  given  cell  size,  the  substitution  of  cadmium  for  zinc 
tends  to  reduce  the  antiferromagnetic  next-nearest-neighbor  interaction,  or  to  increase  the  ferro¬ 
magnetic  nearest-neighbor  interaction,  or  both.  This  conclusion  is  consistent  with  the  relatively 
low  N6el  temperature  in  CdCr2C>4  and  with  the  advent  of  ferromagnetism  in  the  cadmium- 
chromium  chalcogenides. 
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Although  it  would  be  attractive  to  assume  that  the  A-sites  are  directly  involved  in  the  next- 
nearest-neighbor  interactions  of  the  type  Cr-anion-A-site-anion-Cr,  our  results  cannot  rule 
out  an  alternative  explanation  wherein  the  A-sites  play  only  an  indirect  role  in  Cd-anion-anion-Cr 
interactions  by  affecting  the  electron  distribution  of  the  anions. 

N.  Menyuk  R.  J.  Arnott 
K.  Dwight  A.  Woldt 


2.  Reduced  Manganese  Moment  in  Manganese  Chromite 


The  following  is  an  abstract  of  a  paper  which  is  being  presented  at  the  11th  Annual  Confer¬ 
ence  on  Magnetism  and  Magnetic  Materials,  San  Francisco,  California: 

55 

"The  net  magnetization,  magnetic  susceptibility,  Mn  nuclear  magnetic  resonance, 
and  neutron  diffraction  properties  of  MnCr^O^  have  been  reinvestigated.  Our  results 
are  in  general  agreement  with  earlier  findings.  The  neutron  diffraction  pattern  of 

MnCr^O.  evidences  appreciable  deviation  from  a  pure  spiral  configuration.  Further- 

^  4  ++ 
more,  the  susceptibility  data  indicate  that  the  Mn  ions  possess  a  reduced  magnetic 

moment  of  about  4.3  rather  than  the  previously  assumed  value  of  5.0  p^.  Given 

the  hypothesis  of  an  admixture  of  quartet  states  with  quenched  orbital  momentum  into 

6  55 

the  S  free-ion  ground  state,  the  reduced  value  of  the  Mn  hyperfine  field  observed 

in  manganese  chromite  is  consistent  with  this  low  value  for  the  manganese  moment. 


"The  interpretation  of  the  NMR  result  in  terms  of  a  ferrimagnetic  spiral  ground- 
state  spin  configuration  remains  essentially  unchanged,  a  manganese  cone  angle  of 
about  61°  being  indicated.  The  use  of  4.3pR  for  the  Mn++  moment  improves  the  agree¬ 
ment  between  the  theoretical  and  observed  neutron  diffraction  patterns,  but  the  best  fit 
is  still  obtained  with  a  much  smaller  value  for  this  cone  angle,  namely  33°.  Within  the 
spiral  model,  the  observed  magnetic  diffraction  intensities  are  incompatible  with  the 
larger  cone  angle;  the  NMR  findings  are  incompatible  with  the  smaller.  Other  dis¬ 
crepancies  are  also  discussed."  ^  £)w^g^  j  Feinleib 

N.  Menyuk  A.  Woldt 


3.  Temperature  Variation  of  Spin  Wave  Dispersion  Relation 

9  2 

In  previous  experiments,  the  long-range  interaction  <p  >  was  obtained  from  the  ratio  of 

the  coefficient  of  the  quartic  term  F  to  the  coefficient  of  the  quadratic  term  D  in  the  spin  wave 

2  4 

dispersion  relation  =  Dk  —  Fk  by  spin  wave  measurements  at  70  GHz  on  permalloy  at  298°, 
77°,  and  4.2  °K.  Similar  measurements  have  now  been  completed  over  a  temperature  range  of 
4.2°  to  300  °K  at  61.8GHz. 

The  measurements  are  of  interest  because  the  dispersion  relation  and  hence  <p  >  are  ob¬ 
tained  for  the  same  temperature  range  over  which  spin  wave  interactions  and  magnetization 

measurements  have  been  found  to  follow  spin  wave  theory.  In  addition,  the  temperature  depend- 

7  /  2 

ence  of  F,  which  can  be  related  to  a  T  1  term  in  the  temperature  dependence  of  the  exchange 
parameter  D,  is  obtained. 

t  Engineering  Department,  Brown  University. 
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The  experimental  method  and  the  procedure  for  analyzing  the  data  were  the  same  as  pre¬ 
viously  described.  The  apparatus  was  modified  to  include  a  triple  dewar  system  which  allowed 
temperature  control  from  4.2°  to  298  °K  by  means  of  an  exchange  gas.  The  iqagnetic  field  was 
measured  by  means  of  an  NMR  probe. 

It  was  found  that  on  a  63%  Ni-37%  Fe  film,  6950  A  thick,  the  coefficients  of  both  the  quadratic 

and  the  quartic  terms  vary  as  T^'  2  up  to  80  °K  and  approximately  as  T^  ^  up  to  298  °K.  This  is 

the  same  behavior  previously  reported*0  for  the  quadratic  term  using  a  different  permalloy  film 

2  5 /2 

at  9  GHz.  These  results  indicate  that  <p  >  varies  as  T  '  over  the  low  temperature  range,  and 
that  the  range  increases  by  10  percent  from  298°  to  4.2  °K. 

R.  Weber 

P.  E.  Tannenwald 

J.  W.  Burke 


4.  On  MarshalPs  General  Spin  Wave  Theory 

1 1 

Marshall  has  given  the  following  beautifully  simple  and  fairly  general  model  for  treating 
a  system  of  interacting  spin  waves.  Suppose  we  assume  that  spin  waves  are  bosons  described 
by  occupation  numbers  n^(  =  0,  1,  2,  .  .  . ,  k  going  over  a  Brillouin  zone),  and  that  in  the  energy 
eigenstates  (or  approximate  energy  eigenstates)  of  the  system  of  interest,  these  n^  are  good 
quantum  numbers.  Then  the  energy  levels  (or  approximate  energy  levels)  of  the  system  are 
functions  of  the  n^;  for  noninteracting  spin  waves  these  would  be  linear  functions  and  for  low- 
lying  levels  (small  n^)  one  expects  the  expansion 


E({n})  =  2  n  6  —  2  nnc  + 
q  q  q  qp  q  p  qp 


(3) 


to  be  valid.  Following  Marshall,  we  neglect  all  other  types  of  excitations  —  this  means  that 
£  q,  c^,  .  .  .  are  constants  (independent  of  temperature  T)  and  the  energy  levels  E  are  completely 
specified  by  these  n^. 

The  spin  wave  excitation  energy  is  defined  as 


€k  5  <E(.  .  .  ,  nR  +  1,  .  .  .  )  -  E(.  .  . ,  nR,  .  .  .  )> 


(4) 


|njf({n})  exp[-/?E({n})]/2{nj  exp  [-pE({n})  ] ,  and  (3  =  l/(£T).  Assuming 

^  ^  temperature  dependence 


where  <f({n})>  =  2 
finally  that  =  Dk2  —  Fk4>  +  .  .  .  as  k  -*•  0,  Marshall  deduced  the  T 


of  D. 


We  have  simply  followed  this  argument  to  find  the  higher-order  terms  occurring  in  D  as 

well  as  the  temperature  dependence  of  F  (this  involves  the  additional  assumption  that  the  terms 
4 

in  €k  of  order  k  are  spherically  symmetric).  We  obtained 


D  =  Dq(1  -  GT5/2  +  LT7/2  + 


F  =  F  (1  -  HT5/2  +...).  (5) 

The  first  interesting  result  is  the  prediction  that  F  should  also  show  a  T^2  dependence.  This 
was  subsequently  found  to  be  the  case  in  permalloy  films  by  Weber,  Tannenwald,  and  Burke 
(Sec.  IV-C-3).  We  also  noticed  that  the  coefficients  G,  L.  .  .  ,  H.  .  .  are,  surprisingly,  not  linearly 
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independent;  in  fact  an  explicit  relation  between  G,  L,  and  H  follows  from  the  above  considera¬ 
tions  and  is  conveniently  summarized  by 


LT 


7/2 


35F 


GT 


5/2 


=  £T 


4D 


o  £(7/2)  2  H 

'  7  G' 


2  115/2) 


(6) 


where  £(5/2)  =  1.341,  £(7/2)  =  1.127  (Riemann  £-functions).  This  will  allow  an  estimate  of  the 
7  /2 

size  of  the  T  '  term  and  hence  will  indicate  if  this  term  can  be  determined  with  the  present 
experimental  data.  T.A.  Kaplan 

H.  E.  Stanley 


5.  Light  Scattering  from  Spin  Fluctuations  in  Ferro  magnets 


Large  cross  sections  have  been  calculated  for  light  scattering  from  long  wavelength  spin 
fluctuations  in  ferromagnets  near  the  Curie  temperature.  The  spin  fluctuations,  which  appear 
as  magnetization  fluctuations,  modulate  the  optical  dielectric  tensor  comnonents  which  ordinarily 
produce  the  large  Faraday  and  Kerr  effects  below  the  Curie  point. 

bcattering  in  the  present  model  is  caused  by  the  interband  optical  dipole  transitions  be¬ 
tween  3d  and  excited  s  bands,  which  are  perturbed  by  a  spin-orbit  interaction  which  fluctuates 
because  of  spin-spin  interaction.  The  results  can  be  readily  extended  to  the  case  of  localized 
states  in  crystals.  Assuming  that  the  spin-orbit  and  electromagnetic  interaction  can  be  treated 
as  first-order  perturbations,  and  the  spin  fluctuations  due  to  exchange  forces  can  be  decoupled 
from  the  spatial  wave  functions  in  the  unperturbed  system,  the  cross  section  cr  for  light  scattering 
in  the  frequency  range  do;  and  solid  angle  range  dft  is 


d2q 

dcodfi 


r2-^  5F 
0  0)  = 
O 


C(q,  (jOq  —  u)) 


(7) 


where  r  is  the  classical  electron  radius,  and  a>  and  o>  are  the  frequencies  of  the  scattered 
o  '  o  M 

and  incident  light,  respectively.  The  electromagnetic  spin-independent  tensor  form  factor,  JF, 
is  given  in  the  dipole  approximation  by 


l-  l  <fd-fs> 


d,  s 
BZ 


A  -  PdsPsd,  •  A0(*VV(r)xP)d,d 


3  2 


2m  c  lEsd+'ft*  wo  +  lFds^  ^Edd'  +  w  ird'd^ 

(complex  conjugate) 


(8) 


where  A  and  Aq  are  the  unit  polarization  vectors  of  the  scattered  and  incident  light,  with  wave 

vectors  K  +  q  and  K,  respectively.  The  Fermi  factors  for  the  d  and  s  bands,  neglecting  the 

exchange  energy  splitting,  are  f,  .  The  spin- spin  correlation  tensor  Fourier  component  of 

a,  s 

wave  vector  q  and  frequency  o>s  is  C(q,  o)g).  By  the  fluctuation-dissipation  theorem, 

C(q,  o>g)  =  N  Im  x  (q,  <^g)  [coth  (fia>g/2kT)  =f  l]/87TfJig  for  $  0 

where  x  is  the  dynamic  magnetic  susceptibility  and  N  is  the  total  number  of  spins  in  the  crystal. 
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At  low  temperatures,  magnon  and  optical  spin  waves  cause  fluctuations  peaked  about  values  of 

cj  given  by  their  respective  dispersion  curves  o>  (q). 
s  s  — 

Estimates  of  this  cross  section  can  be  obtained  from  the  dielectric  tensor  e  below  the  Curie 
temperature,  namely, 

e(K,  u>)  =  AqI  -  AjI  X  Tr  p°  S  .  (9) 


Then,  with  V  denoting  the  unit  cell  volume,  we  obtain 

3 

i> 

C  (q,  cj  =o>  —  a>) 

zz  —  s  o 


,2 

a*  |  a  y  I 

dcodfi  ,A1  1  A 


2  “0“ 


C  (4tt)' 


(10) 


By  integrating  over  angle  and  frequency  and  using  the  Kramers -Kronig  theorem,  the  total  cross 
section,  in  the  high  temperature  limit,  becomes 


I Alv!2  “o 

*  “  — - 4  N  Xn<T> 


(11) 


where 

X_(T)  =  Re  x(q.  0)  (kT/4fi|)  ~ 


The  Faraday  and  Kerr  effects  on  transmitted  or  reflected  light  are  proportional  to  A^.  For  ex¬ 
ample,  it  is  known  that  for  nickel  and  iron.  Re  A^  “  Im  A^  =*  0.1.  For  (T  —  T^)/T  =*  10~4  (corre¬ 
sponding  to  temperature  control  of  within  0.1  °K),  y  =*  4/ 3  (corresponding  to  measurements  of  the 

2  -4 

long  wavelength  susceptibility),  and  (aq)  =*10  (corresponding  to  an  exchange  interaction  length 

o  —  -3  2  /  2 

a  of  about  10  A),  Eq.  (11)  yields  a  cross  section  cr  =*  10  cm  .  Thus,  an  incident  power  of  1  W/cm 

results  in  a  total  of  a  milliwatt  of  scattered  power  from  a  cubic  centimeter  of  illuminated  crys¬ 
tal.  Scattering  of  a  transmitted  beam  could  be  observed  in  transparent  ferromagnetic  crystals, 
such  as  YIG,  CrBr^,  and  EuS,  whereas  scattering  of  a  reflected  beam  could  be  observed  from 
ferromagnetic  conductors,  such  as  nickel  and  iron. 

The  frequency  spectrum  of  the  spin  fluctuations  is  expected  to  become  sharply  peaked  at 
zero  frequency,  like  Rayleigh  scattering  in  critical  liquids,  because  of  the  large  susceptibility 

near  the  Curie  temperature.  According  to  neutron  scattering  measurements  on  iron,  the  spectral 

2  —3  2 

width  A co  varies  as  D  [q  +  const  (T  —  Tc)].  Since  D  =*  2  X  10  cm  /sec,  the  width  may  be  a  few 
megacycles  or  less,  which  is  suitable  for  laser  "self-beating"  experiments. 

M.  M.  Litvak 
H.  J.  Zeiger 


x  /T  -  T  v 

=' hr1 
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6.  Paramagnetic  Spin  Waves  in  Superconductors 

A  gas  of  interacting  electrons  has  certain  modes  of  oscillation  of  high  phase  velocity  called 
"zero  sound."  The  term  was  coined  by  Landau  who  first  discovered  them.  Since  Landau  proposed 
their  existence,  they  have  remained  experimentally  undetected.  Our  interest  is  in  a  subclass  of 
these  oscillations  which  seems  to  offer  the  greatest  possibility  of  being  detected  in  actual  materi¬ 
als.  This  class  of  modes  is  characterized  by  propagating  spin  density  variations.  We  find  that 
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materials  which  exhibit  superconductivity  are  far  more  likely  candidates  for  the  detection  of 
these  modes;  we  propose  two  experiments,  the  first  modeled  on  a  conventional  thin-film  spin 
wave  resonance  experiment,  and  the  second,  a  more  elaborate  film  transmission  experiment. 

The  second  is  designed  to  observe  an  effect  of  the  electronic  orbital  motion,  a  magnetic  field 
induced  rotation  of  the  wave  fronts  of  the  propagating  spin  wave. 

Because  certain  features  of  the  problem  may  be  solved  without  recourse  to  approximation, 
the  study  of  these  modes  illuminates  some  features  of  collective  mode  theory  in  general.  In 
particular,  while  the  excited  states  of  a  system  described  by  an  Hermitian  Hamiltonian  have 
infinite  lifetimes,  approximate  theory  and  experiments  generally  find  finite  lifetimes  and  nonzero 
resonance  linewidths.  The  difference  between  these  points  of  view  and  the  resolution  of  the  ap¬ 
parent  paradox  may  be  very  clearly  seen  in  studying  the  paramagnetic  spin  wave. 

L.  L.  Van  Zandt 


D.  QUANTUM  TRANSPORT  THEORY 
1.  Pseudo  momentum  of  Phonons 

A  simple  and  general  treatment  of  the  role  of  ftq,  where  q  is  the  wave  vector  of  a  normal 
mode  of  vibration  of  a  lattice,  as  a  pseudomomentum  of  the  phonon  has  been  found  by  consider¬ 
ing  the  transformation  properties  of  the  vibrating  system  under  the  unitary  operator 

C(x)  =  exp  [—  iK  •  x/ft]  ,  (12) 

corresponding  to  the  Hermitian  operator 

K  =  Yj  b-!  b_  ,  (13) 

qK  q* 

where  b  and  b+  are  the  destruction  and  creation  operators  for  the  phonons. 

For  the  anharmonic  lattice,  i.e.,  interacting  phonons,  it  is  easy  to  see  that  the  Hamiltonian 

is  invariant  under  the  unitary  transformation  C(am>,  where  am  is  a  lattice  vector,  i.e., 

C(a  )  HC_1(a  )  =  H.  Thus,  C(a  )  is  a  constant  of  the  motion,  [C(a  ),  H]  =  0.  If  now  we  write 
mm  m  m 

H  =  Hq  +  H^,  where 


h  =  y  (b*  t>_  +  4) 

qK  qK  qK 


(14) 


qK 


is  the  Hamiltonian  for  the  harmonic  part  of  the  vibrations  and  H^  the  phonon-phonon  interaction, 
we  have 


[C(a  ),  H,]  =  0 

i  m" 


(15) 


since  C(am)  clearly  commutes  with  HQ. 
Eq.  (15)  gives 


In  a  representation  in  which  all  b_^  b_^ 


qK  qK 


are  diagonal. 
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exp 


—  4-  K  {n  }  •  a 

fi  —  m 

q* 


]-exp[—  L  kR}  •  am] 


<{n_  }  |H1 1  {n^  }>  =  0  .  (16) 

q  k  q*r 


Thus,  for  all  processes  possible  due  to  the  phonon-phonon  interaction  H^,  we  must  have  the 
selection  rule 

K  {n_  }  =  K  {n^  }  +  ft  b_  ,  (17) 

q  k  q/c  n 


where  b_^  is  a  reciprocal  lattice  vector.  Thus,  we  have 
n 


X  3n-.  =  X  +  ,  (18) 

QK  qk  n 

qK  M  qK  H 

which  exhibits  the  pseudomomentum  character  offiq  for  each  phonon. 

A  similar  treatment  can  be  given  for  the  harmonic  crystal  interacting  with  an  external  sys¬ 
tem  which  cannot  distinguish  between  the  crystal  particles.  This  system  can  be  proved  to  be 

invariant  under  C(a  )  and  a  simultaneous  translation  of  the  external  system  by  a  .  This  leads 
m  J  J  m 

in  an  analogous  way  to 


p+  2j^qn—  =  P 1  +  +  ft  b_^ 

qK  —  qK  n 


q* 


qK 


(19) 


where  p  and  pT  are  the  momenta  of  the  external  system  in  its  "initial"  and  "final"  states. 

Note  that  the  transformation  by  C(am>  does  not  translate  the  particles.  The  operator  that 
generates  the  translation  is,  of  course,  the  total  momentum  of  the  lattice.  This  leads  to  the 
familiar  conservation  of  the  total  momentum,  in  addition  to  Eqs.  (18)  and  (19).  The  main  point 
here  has  been  anticipated  in  the  work  of  G.  Sussmann.^ 

P.  N.  Argyres 


2.  Polaron  Energy  Spectrum 

The  energy  spectrum  of  a  polaron  —  an  electron  interacting  with  the  longitudinal  optical 
phonons  of  a  polar  crystal  —  has  been  studied,  with  particular  attention  given  to  the  region  where 
the  polaron  excitation  energy  closely  approaches  the  energy  of  a  longitudinal  optical  phonon.  In 
this  region,  existing  theories  of  the  polaron  spectrum  are  inadequate;  in  particular,  the  usual 
Rayleigh-Schrodinger  perturbation  theory  is  shown  to  be  inconsistent. 

A  self-consistent  weak  coupling  theory  and  a  variational  theory  of  the  polaron  spectrum  which, 

for  small  coupling,  reduces  to  this  weak  coupling  theory,  have  been  developed. 

13  14 

Schultz  and  Whitfield  and  Puff  have  conjectured  that  the  polaron  energy  E(p)  bends  over 
and  becomes  horizontal  as  the  polaron  momentum,  p,  approaches  from  below  the  value  at  which 
the  polaron  excitation  energy,  E(p)  —  E(0),  becomes  equal  to  the  optical  phonon  energy.  Using 
the  above-mentioned  weak  coupling  method,  this  conjecture  has  been  verified  to  lowest  order 
and  next  higher  order  in  the  coupling  constant.  n  m  T 


3.  Quantum  Theory  of  Kinetic  Equations  for  Electrons  in  Random  Impurities 

We  have  considered  the  problem  of  deriving  kinetic  equations  for  electrons  under  the  action 
of  arbitrary  space-time  varying  driving  fields  and  randomly  distributed  impurities  from  a  fresh 
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point  of  view.  We  first  derive  a  kinetic  equation  for  the  whole  one-electron  density  operator  by 
eliminating  the  dependence  on  the  positions  of  the  impurities,  using  a  general  procedure  reported 
earlier  by  Argyres.1^ 

From  this  general  kinetic  equation  we  can  then  easily  derive  special  kinetic  equations  for 
the  subsets  of  the  one-electron  density  matrix  appropriate  to  special  situations.  For  example, 
the  usual  Boltzmann  transport  equation  for  the  distribution  function  is  derived,  without  the  ad¬ 
ditional  assumptions  or  theorems  necessary  in  earlier  theories.  The  appropriate  kinetic  equa¬ 
tion  when  a  quantizing  magnetic  field  is  present  is  also  easily  deduced. 

The  theory  of  power  absorption  by  direct  and  indirect  processes  is  developed  on  the  basis 
of  the  general  kinetic  equation. 

The  case  of  an  inhomogeneous  electromagnetic  field  is  considered  in  detail,  with  emphasis 
on  the  requirements  of  gauge  invariance  and  the  conditions  necessary  for  the  commonly  assumed 
form  of  the  Boltzmann  equation  for  the  space-varying  distribution  function. 

A  full  report  of  the  work  is  being  prepared  for  publication. 

P.  N.  Argyres 
E.  S.  Kirkpatrick 

4.  Quasi-Photons  in  Charged  Media 

The  interaction  of  photons  with  charged  particles  leads  in  the  zeroth  approximation  to  quasi¬ 
photons  with  dispersion  given  by 

2  2  2  2 

o>k  =  c  k  +  Wp  ,  (20) 

where  o)p  is  the  plasma  frequency  of  the  charged  particles.  This  relation  has  recently  been 
derived by  the  method  of  Bogoliubov  and  Tyablikov  and  also  by  the  method  of  the  equation  of 
motion. 

We  have  proved  relation  (20)  in  a  simple  way  without  any  transformation,  by  first  writing 
the  Hamiltonian  for  the  system  not  in  terms  of  the  creation  and  destruction  operators  for  the 
bare  photons  and  the  charged  particles,  but,  instead,  in  terms  of  the  complex  normal  coordinates 
of  the  fields.  It  is  then  immediately  obvious  that  in  addition  to  the  "potential"  energy  of  the  bare 
photons 

£  c2k2A  +  A  ,  (21) 

kA  k*  k* 

there  is  a  similar  term  in  the  Hamiltonian  due  to  the  A2-term  of  the  interaction  of  the  radiation 
field  with  the  mean  number  of  the  charged  particles.  This  term  can  be  written  in  terms  of  the 
complex  normal  coordinates  as 


Z 

kX 


2  A  +  A 
Ci)  A  A 

p  kX  kX 


(22) 


which  combines  with  the  previous  term  of  the  bare  photons  to  give  the  renormalized  photons 
with  dispersion  given  by  Eq.  (20).  p  N  Argyres 
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